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ABSTRACT 

We have conducted a large-scale survey of 120 planetary nebulae (PNe) to 
search for the near-infrared emission lines [Kr III] 2.199 and [Se IV] 2.287 /im. 
The neutron(n)-capture elements Se and Kr may be enriched in a PN if its 
progenitor star experienced s-process nucleosynthesis and third dredge-up. In 
order to determine Se and Kr abundances, we have added these elements to the 
atomic databases of the photoionization codes Cloudy and XSTAR, which we use 
to derive ionization correction factors (ICFs) to account for the abundances of 
unobserved Se and Kr ions. However, much of the atomic data governing the ion- 
ization balance of these two elements are unknown, and have been approximated 
from general principles. We find that uncertainties in the atomic data can lead to 
errors approaching 0.3 dex in the derived Se abundances, and up to 0.2-0.25 dex 
for Kr. To reduce the uncertainties in the Kr ionization balance stemming from 
the approximate atomic data, we have modeled ten bright PNe in our sample, 
selected because they exhibit emission lines from multiple Kr ions in their op- 
tical and near-infrared spectra. We have empirically adjusted the uncertain Kr 
atomic data until the observed line intensities of the various Kr ions are ade- 
quately reproduced by our models. Using the adjusted Kr atomic data, we have 
computed a grid of models over a wide range of physical parameters (central star 
temperature, nebular density, and ionization parameter), and derived formulae 
that can be used to compute Se and Kr ICFs. In the second paper of this series, 
we will apply these ICFs to our full sample of 120 PNe, which comprises the first 
large-scale survey of n-capture elements in PNe. 
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stars: AGB and post-AGB — atomic data — infrared: general 

1. INTRODUCTION 

The compositions of planetary nebulae (PNe), the ejected envelopes of evolved low- 
and intermediate-mass stars (M = 1-8 Mq), are affected by past nucleosynthesis in their 
progenitor stars. In particular, these objects are important sources of C, N, and roughly half 
of the neutron(n)-capture element (atomic number Z > 30) nuclei in the Universe (Busso et 
al. 1999, hereafter BGW99). 

Carbon and n-capture elements are produced during the thermally-pulsing asymptotic 
giant branch (AGB) phase. Trans-iron nuclides are formed by slow ra-capture nucleosynthesis 
(the "s-process" ) , in which iron peak "seed" nuclei in the region between the H- and He- 
burning shells capture free neutrons released primarily by the reaction ^^C{a,n)^^0 (or, in 
more massive AGB stars, ^^Ne(«,n)^^Mg). The seed nuclei undergo a series of n-captures 
and /3-decays that transform them into heavier elements (Kappeler et al. 1989; BGW99; 
Goriely & Mowlavi 2000; Lugaro et al. 2003; Herwig 2005). 

The intershell material, enriched in C (from partial He burning) and n-capture elements, 
is conveyed to the stellar surface by convective mixing, a process denoted "third dredge-up" 
(TDU) to differentiate it from earlier dredge- up events during the red giant branch and early 
AGB phases. TDU is believed to occur in stars with initial masses greater than 1.2-1.5 M©, 
and is characterized by recurrent nuclear runaways of the He-burning shell (i.e., thermal 
pulses) separated by periods of convective mixing and transport to the surface (Iben & 
Renzini 1983; BGW99; Herwig 2005; Straniero et al. 2006). The enriched material is then 
released into the surrounding interstellar medium (ISM) by stellar winds and ultimately PN 
ejection. 

Abundance determinations in PNe can be used to investigate the occurrence of the 
s-process and TDU in PN progenitor stars, and their importance to the Galactic chemical 
evolution of n-capture elements. However, the low abundances of n-capture elements (< 10~^ 
relative to H; Asplund et al. 2005) cause their emission lines to be weak, and these elements 
were not seriously considered to be detectable in PNe until just over a decade ago. 

It was not until 1994 that n-capture element lines were identified in the spectrum of a PN 
(Pequignot & Baluteau 1994, hereafter PB94). These authors detected emission lines of Kr 
{Z = 36), Xe {Z = 54), and possibly other n-capture elements in a deep optical spectrum 
of the bright PN NGC 7027. Using approximations to the unknown Kr and Xe collision 
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strengths (which have since been calculated by Schoning 1997 and Schoning &: Butler 1998) 
and corrections for the abundances of unobserved ionization stages, these authors found that 
Kr and Xe are highly enriched in NGC 7027, providing evidence for the occurrence of the 
s-process and TDU in its progenitor star. 

The study of PB94 led Dinerstein (2001) to realize that two long- unidentified emission 
lines seen in the K band spectra of several PNe are in fact fine-structure transitions of 
[Kr III] and [Se IV]. She found that the observed Kr and Se {Z = 34) line strengths in 
IC 5117 and NGC 7027 are consistent with enhanced abundances, and furthermore proposed 
that the presence of these lines in some PNe but not others indicates a spread in s-process 
enrichments among Galactic PNe. 

Until now, the only other studies of n-capture element abundances in PNe were those 
of Sterling et al. (2002) and Sterling & Dinerstein (2003a), who detected Ge {Z = 32) 
in absorption against the UV central star continua of six PNe with the Far Ultraviolet 
Spectroscopic Explorer (FUSE). They derived Ge abundances in five of the objects, and 
found it to be enhanced in four PNe by factors of > 3-10, depending on the amount of 
Ge depletion into dust. For one of these objects, SwSt 1, Sterhng et al. (2005) determined 
the fine of sight Fe abundance from the FUSE spectrum, and found it to be only slightly 
depleted ([Fe/S] = —0.35 ±0.12). This implies that Ge is also negligibly depleted in SwSt 1, 
and that [Ge/S] = 0.72 ± 0.06, a factor of five enrichment relative to solar. 

In all, these studies resulted in n-capture element abundance determinations in only 
seven PNe. Such a small sample of objects divulges little information regarding the s- 
process and dredge-up histories of PN progenitors as a population, or their role in the 
Galactic chemical evolution of trans-iron species. 

Therefore, we have undertaken a survey of Galactic PNe to search for the two near- 
infrared (NIR) lines identified by Dinerstein (2001), [Kr III] 2.199 and [Se IV] 2.287 /xm. Se 
and Kr are useful tracers of s-process enrichments and TDU in PN progenitor stars for two 
reasons: (1) Neither is expected to be depleted into dust, as Kr is a noble gas and Se has 
not been found to be significantly depleted in the diffuse ISM (Cardelli et al. 1993); and (2) 
the detection of the NIR lines of these elements in several PNe (e.g., Geballe et al. 1991; 
Hora et al. 1999; Lumsden et al. 2001) indicates that they may be observable in a significant 
fraction of Galactic PNe. Furthermore, theoretical models predict that Se and Kr can be 
significantly enriched by s-process nucleosynthesis in AGE stars (e.g., Gallino et al. 1998; 
Goriely & Mowlavi 2000; Busso et al. 2001). Therefore, if a star experiences the s-process 
and efficient TDU, Se and Kr will be enriched in the PN it produces. 

We have observed 103 PNe in the K band, and utilize hterature data to expand our 
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sample to 120 objects. This is the first large-scale survey of n-capture elements in PNe, and 
comprises the first broad characterization of s-process enrichments in PNe as a population. 
We have detected Se and/or Kr in 81 of the 120 PNe in our sample, which increases the 
number of PNe with known n-capture clement abundances by more than a factor of ten. 
Preliminary results from our survey have been presented in Sterling & Dinerstein (2003b; 
2004; 2005a,b; 2006) and Sterling (2006a). 

In order to derive elemental Sc and Kr abundances from our data, it is necessary to 
correct for the abundances of unobserved ionization stages. We have added these elements 
to the photoionization codes Cloudy (Ferland et al. 1998) and XSTAR (Kallman & Bautista 
2001; Bautista & Kallman 2001), which we use to derive widely applicable formulae for 
calculating Se and Kr ionization correction factors (ICFs). 

This task is complicated by the general lack of atomic data governing the ionization 
balance of Se and Kr (e.g., photoionization cross-sections and rate coefficients for various 
recombination processes). In fact, these data are unknown for almost all n-capture elements. 
Therefore, we have used several approximations to calculate the atomic data of Se and Kr. 
We have modeled ten PNe in our sample, which were selected because they exhibit emission 
lines from multiple Kr ions in their optical and NIR spectra. This enables us to adjust 
the uncertain Kr atomic data until the observed line intensities of the various Kr ions arc 
satisfactorily reproduced by our models. Unfortunately, no other Se ions have been detected 
in PNe, and thus we are unable to empirically correct the uncertain Se atomic data. With 
the adjusted Kr atomic data, wc have run a grid of Cloudy models spanning a wide range of 
physical parameters (central star temperature, nebular density, and ionization parameter). 
We extracted the fractional abundances of Se and Kr ions from this grid, and searched for 
correlations between these and the ionic fractions of more abundant elements. The best 
correlations have been fit with analytical functions, which can be used to compute Se and 
Kr ICFs. This paper, the first of two, describes these model calculations and presents the 
best analytical fits to the Se and Kr ICFs. In the second paper of this scries (Sterling & 
Dinerstein 2007, hereafter Paper II), we will use our model-derived ICFs to determine Se 
and Kr abundances for our full sample of objects. 

The outline of this paper is as follows: in §2, we briefly describe our observations 
and supporting data from the literature for the 10 objects we model in detail. In §3, we 
present our models of these objects, empirically adjust the uncertain Kr atomic data, and 
discuss how uncertainties in the adopted atomic data affect the derived abundances. In 
§4, we present the grid of Cloudy models from which we derive formulae to correct for 
unobserved ionization stages of Se and Kr. We discuss our results and provide concluding 
remarks in §5. Two appendices are also included. Since the optical and NIR fines of Se 
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and Kr have been frequently misidentified in the hterature, in Appendix A we investigate 
alternative identifications to the lines we consider in this study. In particular, we remove 
important contaminants to the Se and Kr lines, and systematically rule out other possible 
identifications. In Appendix B, we describe the Se and Kr atomic data we have adopted in 
our photoionization models. 

2. OBSERVATIONS AND LITERATURE DATA 

2.1. Infrared Spectra 

We have observed 103 PNe in the K band with the CoolSpec spectrometer (Lester et 
al. 2000) on the 2.7-m Harlan J. Smith telescope at McDonald Observatory. In addition, 
we utilize K band PN spectra from the literature to expand our sample to 120 objects. In 
Paper II we will fully describe the observational data and reductions, and here we present 
only a brief synopsis. 

Each PN was observed from 2.14-2.30 /im with a 2'.'7 x 90" slit and 75 1/mm grating, at a 
resolution of R ~500 (FWHM = 0.004 /xm at 2.20 /im). All data were reduced using 
For each two-dimensional spectrum, we performed dark subtraction, cosmic ray removal, 
and flat-fielding before extracting the one-dimensional spectrum. The ID spectra were then 
wavelength calibrated and dispersion-corrected using Ar lamp spectra. For each PN, we 
observed at least one AO standard star at a similar airmass, which we used to response- 
correct the spectrum, remove telluric features, and perform flux calibrations. 

[Kr III] 2.199 and [Se IV] 2.287 /im are resolved from nearby features of other species, 
except in PNe exhibiting vibrationally-excited molecular hydrogen. In these objects, the 
Kr and Se lines may be contaminated by H2 3-2 S'(3) 2.201 and H2 3-2 S{2) 2.287 /im, 
respectively. As shown in Appendix A. 4, these lines are not important contaminants to the 
[Kr III] and [Se IV] fluxes of the ten PNe considered in this paper. 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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2.2. Optical Spectra and Source Selection 

We have selected 10 PNe from our sample that exhibit emission from [Kr IV] (and in 
one case, [Kr V]) in their optical spectra, in addition to optical or NIR [Kr III] emission lines. 
The additional observed ionization stages reduce the uncertainties in correcting for unseen 
Kr ions, and enable us to empirically adjust the approximate Kr atomic data with the aid 
of photoionization models. Unfortunately, no other Se ions have been clearly identified in 
PNe, and hence we are unable to optimize the Se atomic data. (A possible exception is 
[Se III] A8854.2, identified by PB94 in NGC 7027. However, this line was undetected in the 
higher quality spectrum of Zhang et al. 2005, and hence we consider PB94's identification 
doubtful.) 

The optical data references for these ten PNe are given in Table 1, along with the spectral 
resolution, extinction coefficients, and the sources of the NIR data and (when applicable) 
other useful information, such as elemental abundances and UV and IR line intensities. The 
NIR line fluxes were dereddened using the extinction coefficients (Table 1) and adopted 
extinction law of the optical references. We then applied an aperture correction to the 
NIR data (to account for the different slit widths and positions of the optical and NIR 
measurements) by forcing the H I Bi'j/B.p intensity ratios into agreement with the theoretical 
ratio at the Tg and Ue determined from the optical data, using Table B.9 of Dopita & 
Sutherland (2003). 

The optical Kr line intensities are given in Table 2, where the cited 1-a errors are the 
fiux uncertainties indicated in the optical data references. Table 2 also fists the dereddened 
and aperture-corrected NIR Se and Kr fine intensities of these objects. The intensities have 
been corrected for contaminating blends as described in Appendix A. 

3. PHOTOIONIZATION MODELS 

In order to derive Se and Kr abundances that are as accurate as possible, we have used 
two independent photoionization codes. Cloudy (Ferland et al. 1998) and XSTAR (Kallman 
& Bautista 2001; Bautista & Kallman 2001), to model the ionization structure of Se and Kr. 
We have expanded the atomic database for each of these codes to include Se and Kr, and 
performed detailed spectral fits of the ten PNe listed in Table 1. 

The Se and Kr atomic data that we have utifized in Cloudy and XSTAR are discussed 
in Appendix B. Much of the atomic data controlling the ionization balance of these elements 
is unknown, in particular the photoionization (PI) cross-sections and rate coefficients for 
radiative recombination (RR), dielectronic recombination (DR), and charge transfer (CT). 
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As explained in Appendix B, we have used various approximations to calculate values for 
the cross-sections and rate coefficients of these atomic processes. 

The atomic data has been entered into Cloudy version C06.01 and XSTAR version 
2.1kn3. Cloudy loops through all elements up to the maximum atomic number in its 
database, and therefore it was necessary to enter atomic data for the remaining elements be- 
tween Zn (the heaviest element in vCOG.Ol) and Kr. We computed PI cross-sections and RR 
rate coefficients for Ga, Ge, As, and Br using the same methods described in Appendix B, 
but assumed the CT, DR, and collisional ionization rate coefficients to be the same as the 
corresponding charge states of Se or Kr. These four elements were turned off in all of our 
Cloudy models, and therefore do not affect our results. In contrast, it was not necessary to 
enter data for these intervening elements in XSTAR. 

Se and Kr are trace elements in PNe, even when enriched by the s-process. As such, 
they are not expected to significantly affect the thermal or ionization structure and predicted 
spectra of nebulae. For each of the models described in this section, we also ran the corre- 
sponding unaltered versions of Cloudy and XSTAR (without Se and Kr) with the same input 
parameters to test this assumption. In all cases, we found negligible differences between the 
models with and without Se and Kr, indicating that our additions to the atomic databases 
did not introduce any unforeseen effects in the functionality of these codes. 

In the following, we present detailed models of ten PNe in which the input parameters 
were optimized to best reproduce the observed UV, optical, and IR spectra. We discuss 
uncertainties in the derived Se and Kr abundances in these nebulae arising from a variety of 
sources, including the atomic data, and how these will affect the derived abundances of our 
Ml sample of 120 PNe. 

3.1. Cloudy Models 

We performed detailed spectral fits to the UV, optical, and IR line intensities of the 
ten PNe in Table 1, in order to reproduce their ionization structure. The main objective of 
this exercise was to compare the predicted Kr line strengths to the observed ones, and to 
make necessary adjustments to the uncertain Kr atomic data (Appendix B) so that these 
lines are best reproduced by our models. Once agreement was achieved, we incorporated the 
optimized Kr atomic data into our grid of Cloudy models (§4), which we use to derive ICFs 
for Se and Kr. 

We made a number of simplifying assumptions in these models. First, spherical geometry 
was assumed, with filling and covering factors (the fraction of material in a solid angle of 47r 



-8- 



steradians around the central star) of unity. Second, we assumed the nebulae have constant 
densities. Since this produced reasonably good fits to the observed spectra, we did not 
attempt to use other density laws. Finally, we do not include dust physics in our models of 
single objects, aside from the depletion of refractory elements into grains. This is likely to 
result in an inaccurate treatment of the thermal balance of the nebulae, due to the influence 
of photoelectric heating by dust (Dopita & Sutherland 2000). For this reason, the exclusion 
of dust may also affect the derived central star temperatures. In §4, we show that ignoring 
dust physics has a neghgible effect on the predicted ionization balance of Se and Kr. 

These simplifications arc not likely to be important for the purposes of our study. Our 
goal is not to produce physically accurate models of these PNe, but rather to establish the 
ionization balance of the nebulae in order to adjust the approximate Kr atomic data. For 
this reason, our model-derived abundances for elements other than Se and Kr should not be 
regarded as improvements over those of the references listed in Table 1. 

For the input parameters of each model, we initially adopted the abundances and den- 
sities reported in the hterature (see Table 1). It is also necessary to specify the stellar flux 
distribution, inner and outer radii of the nebular material, and ionizing luminosity of the 
central star. 

In all cases, we utihzed NLTE stellar atmospheres including metal fine blanketing, com- 
puted for hot compact stars such as the central stars of PNe (Ranch 2003). Solar abundances 
were assumed for the model atmospheres, since the modeled PNe are all Galactic disk objects 
hkely to have near-solar abundances. Since the Rauch (2003) grid only extends from 50,000 K 
to 190,000 K, we were unable to use it to model IC 418, which has a very cool central star. 
T. Rauch (2005, private communication) kindly computed models with Tcs = 30, 000 and 
40,000 K and log(5') = 3.4, which we use to model IC 418. The Rauch stellar atmospheres 
do not include mass-loss, which can significantly alter the ionizing fiux of a PN central star 
(Kudritzki & Pols 2000). We have tested the dependence of our results on this factor by 
modehng IC 418 with a stellar atmosphere including mass-loss (Koesterke et al. 2007, in 
preparation), and find that incorporating mass- loss does not significantly affect the derived 
nebular abundances (§3.1.1). 

For the nebular radii and stellar luminosities, we used as initial parameters the quantities 
derived from the photoionization models of the Hyung and AUer references cited in Table 1, 
in addition to Hyung et al. (1994b) for IC 418, Hyung & AUer (1997a) for NGC 6741, Hyung 
et al. (1997) for NGC 6884, Keyes et al. (1990) for NGC 7027, and Hyung & AUer (1997b) 
for NGC 7662. We also adopted the stellar gravities of these references for the Rauch NLTE 
models, except for IC 418, where we used log{g) = 3.4 rather than the value of 2.9 computed 
by Hyung et al. (note that Mendez et al. 1992 found log{g) = 3.45 for this star). In the 
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case of NGC 6826, we used the nebular radius, Tes, and stellar luminosity of Mendez et al. 
(1992) as initial parameters. Instead of log(5') = 4.0, as computed by those authors, we used 
log(gi) = 5.0, which is the minimum gravity in the grid of Ranch (2003) NLTE atmospheres. 

With initial values for all nebular parameters in hand, we varied the abundances, T^s, 
riB., and in many cases the outer radii and stellar luminosity until a good fit between the 
predicted and observed spectra was obtained. We used a mixture of the subplex optimization 
routine in Cloudy and manual manipulation of the parameters to achieve these fits. The 
input parameters for all our models are given in Tables 3 and 4 (with the solar composition 
in the last column for comparison). In Tables 5, 6, and 7, wc display the predicted versus 
observed line intensities (relative to H/3) for lines used in the spectral fits. The intensities 
are all on the scale /(H/5) = 100. For most lines, the predicted intensities are within 25% of 
the observed ones. In some cases, there are larger discrepancies; these are hkely due to the 
simplifying assumptions discussed above. 

With the Kr atomic data we originally calculated, we found that the predicted [Kr III] 
fines were too strong, and the [Kr IV] fines too weak relative to the observed intensities. In 
addition, the [Kr V] A6256.1 intensity was underproduced in NGC 7027, the only PN of our 
sample in which this line was clearly detected. Increasing the Kr++ PI cross-section derived 
from Equation (Bl) by a factor of two, and decreasing that of Kr^"*" by a factor of four leads 
to better agreement in the relative line intensities of [Kr III] , [Kr IV] , and [Kr V] . 

We adopt these optimized PI cross-sections (listed in Table B.2) in our models. The 
predicted Kr fine strengths in Tables 5-7 are from models using the adjusted Kr"^"*" and 
Kr^"*" PI cross-sections. As seen in these tables, our models do not predict the Kr ionization 
balance to be systematically too high or low. The differences between the predicted and 
observed Kr line intensities can be attributed to flux uncertainties (~30-40%) for all but 
two objects. In NGC 6884, [Kr III] A6826.7 is predicted to be too strong by 60%, but as 
discussed in Appendix A. 3, the flux of this line is uncertain due to possible blending witfi 
He I and telluric OH features. In IC 418, [Kr III] 2.199 fim is predicted to be too strong, 
while [Kr IV] fines are too weak. We were unable to reproduce the observed [Kr IV] line 
intensities in this object except with a higher T^s, which led to a poor fit to the rest of the 
spectrum. 

Note that the derived Se and Kr abundances given in Tables 3-4 are often higher than 
the solar Se and Kr abundances (3.36±0.04 and 3.28±0.08, respectively) of Asplund et al. 
(2005), particularly in the case of Kr. In many cases, these enrichments are too large to 
be explained by the observed star-to-star scatter in the initial abundances of n-capture 
elements in stars of near-solar metallicity (Travaglio et al. 2004; Simmerer et al. 2004). 
Such large abundance enhancements must instead be due to self-enrichment by 5-process 
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nucleosynthesis (BGW99). In addition, theoretical s-process models predict Kr to be more 
enriched than Se by 0.0-0.5 dex (Gallino et al. 1998; Goriely & Mowlavi 2000; Busso et al. 
2001), as we have found for most of these PNe. We believe that the subsolar Se abundances 
we derive for NGC 6790 and NGC 6826 are likely to be reflections of the overall lower 
metallicities of these objects rather than any type of Se depletion (e.g., into dust), based on 
the subsolar O, Ne, and Ar abundances determined for these objects by Liu et al. (2004b). 
We reserve a full discussion of the s-process enrichments in these objects for Paper II. 



3.1.1. Dependence of IC 418 Model on Stellar Mass-Loss 

In order to test the sensitivity of our results to the effects of stellar winds on the 
ionizing fluxes of the central stars (Kudritzki & Pols 2000), we have modeled IC 418 with 
a stellar atmosphere including mass-loss, computed for us by L. Koesterke (2005, private 
communication; see also Koesterke et al. 2007, in preparation). This stellar wind model was 
computed by fitting a synthetic spectrum to International Ultraviolet Explorer (lUE) low- 
dispersion SWP spectrum retrieved from the Multi-Mission Archive at the Space Telescope 
Science Institute (MAST). The low-dispersion SWP datasets taken with the large aperture 
(SWP03177, SWP03178, SWP03810, SWP06651, SWP43170, SWP43171, SWP46377, and 
SWP46378), were weighted by exposure time and co- added to increase the signal-to-noise 
ratio. A Ranch NLTE atmosphere with no mass-loss was used initially, and the mass- 
loss rate M was increased until the resonance line P Cygni profiles in the lUE spectrum 
were reproduced. Two stellar models were computed, with T^s = 35,000 and 40,000 K, 
log(5') = 3.4, and M — 10~^ Mq yr~^, which we interpolated to produce a Cloudy model of 
IC 418. 

Our best-fit values for the stellar and nebular properties using the Koesterke stellar 
atmosphere are quite similar to those with the Rauch (2003) atmosphere (see Tables 3 and 
5), with the T^s agreeing to within 5% and the nebular abundances to within 30% with 
the exception of O, Ne, and Kr. The reason for these discrepant abundances lies in the 
flux distributions of the two stellar atmospheres. The mass-losing atmosphere of Koesterke 
emits more O"*" ionizing photons, and therefore produces O"'""'" more easily than the Rauch 
(2003) model. This results in a better fit to the observed O lines, and it is not necessary to 
increase the O abundance to produce agreement with the observed [O III] line intensities, 
as is necessary with the Rauch atmosphere. Similarly, the wind atmosphere ionizes Kr^"*" 
more easily than the Rauch atmosphere, allowing the [Kr IV] lines to be fit with a lower Kr 
abundance. In contrast, the Koesterke atmosphere does not produce enough Ne"*" ionizing 
photons, which does not permit a good fit to the Ne line intensities. It was necessary to 
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increase the Ne abundance until the of the [Ne II] and [Ne III] hne fits was minimized 
(this does not significantly affect the Kr ionization balance) . 

The only other modeled PNe which exhibit wind signatures in their lUE spectra are 
NGC 6572 and NGC 6826. Given the general agreement between Cloudy models of IC 418 
using stellar atmospheres with and without mass-loss, we have not attempted to model these 
other objects with atmospheres including stellar winds. 

3.1.2. Error Analysis 

We have estimated uncertainties for the derived Teflf and the Se and Kr abundances of 
each modeled PN. We do not present formal error bars for the abundances of other elements, 
since it was not our objective to determine accurate abundances for those species. 

Our error estimates for Se and Kr include uncertainties in the line intensities, Teff, 
and hydrogen density nn- First, we estimated the uncertainties in the adopted Teflf by 
manually adjusting this parameter in the models until the nebular temperature diagnostics 
and/or ionization balance were no longer acceptably reproduced (to within 50-100%). These 
uncertainties are given in Tables 3-4. We then re-derived the Se and Kr abundances with 
the maximum and minimum allowed Teg, to determine their sensitivity to uncertainties in 
the assumed T^s- Next, we computed the Se and Kr abundances after changing the density 
by a factor of two (also adjusting the outer radii of the nebulae so that the models did not 
become too optically thin or thick) . Finally, we estimated the errors due to uncertainties in 
the Se and Kr line intensities. For Kr, this involved adjusting the abundance until the most 
discrepant line intensities in Tables 5-7 were reproduced; if all Kr lines were predicted with 
good accuracy (e.g., NGC 7027), we simply assumed errors of 30% in the line strengths. On 
the other hand, only one line of Se was observed for each object. The flux errors for this 
line (Table 2) lead to uncertainties of only 5-25% in the Se abundance, likely resulting in 
underestimated formal errors to the Se abundances. 

In Table 8, we compare the Se and Kr abundance errors (in dex) from each source of 
uncertainty we considered. In general, the line strengths are the greatest source of uncer- 
tainty for the Kr abundances. However, the errors from T^q uncertainties are large in the low 
excitation PNe IC 418 and NGC 6826, where the Kr ionization balance is strongly sensitive 
to the stellar temperature. For Se, the abundance errors have comparable contributions from 
uncertainties in the line strengths, Tes, and nn. We added the errors in quadrature to obtain 
the formal (l-cx) abundance uncertainties given in Tables 3, 4, and 8. 
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3.1.3. Abundance Errors From Atomic Data Uncertainties 

One of the most important sources of uncertainty in our derived Se and Kr abundances is 
the approximate atomic data we have adopted. In order to gain insight into the magnitude of 
the abundance uncertainties introduced by the use of the approximate atomic data, we have 
run Monte Carlo simulations for four of our Cloudy models, in which the threshold PI cross- 
sections and rate coefficients of RR, CT, and low-T DR have been varied separately. The 
Monte Carlo simulations were run by using a Gaussian random number generator included 
in the Cloudy distribution. 

As described in Appendix B, we have used the Kramers (1923) formula to derive thresh- 
old PI cross-sections for Se and Kr ions. Gould (1978) showed that this approximation is 
accurate to within a factor of two for first and second row elements of the Periodic Table. 
Since Se and Kr are significantly more complex systems, the Kramers formula may not be 
as accurate. Therefore, we have varied the Se and Kr PI cross-sections in Cloudy by allow- 
ing each threshold cross-section to randomly take values within a Gaussian distribution of 
FWHM 0.5 dex (a factor of three uncertainty, 1-a) around the value given in Table B2. We 
did not attempt to investigate deviations from the adopted shape of the cross-section as a 
function of energy (a power law of index —2), and this remains an important uncertainty. 
We varied the rate coefficients for CT, low-T DR, and RR by the same amount. While the 
RR rate coefficients were derived in a self-consistent manner using the Milne relation and 
the PI cross-sections, the CT and low-T DR rate coefficients could not be estimated from 
first principles. Instead, we assumed that Se and Kr ions have the same rate coefficients for 
these processes as similar charge states of nearby elements. The 0.5 dex 1-cr uncertainties 
we adopt for the CT and low-T DR rate coefficients are based on the expected uncertainties 
of Landau-Zener CT calculations (Kingdon & Ferland 1996) and the error bars of the low-T 
DR rate coefficient estimates in Cloudy (Ferland et al. 1998), but may underestimate the 
actual uncertainties for Se and Kr ions. 

We chose four PNe to illustrate the errors from atomic data uncertainties: IC 418, 
NGC 6741, NGC 6884, and NGC 7027. These objects were selected because they span a 
large range of nebular density and excitation. For each of these PNe we ran 100 Cloudy 
models, in which the PI cross-sections and rate coefficients for RR, CT, and low-T DR were 
separately varied for all Se and Kr ions. The atomic data of each ion was varied independently 
of other ions. 

In Table 9, we display the results of our Monte Carlo simulations, showing the errors 
(in dex) of predicted line intensities for each observed ion. The intensity errors refiect 
uncertainties in the abundance of the parent ion, which in turn affects the derived elemental 
abundances. It can be seen that errors in the PI cross-sections cause the greatest uncertainties 
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in the predicted line intensities, followed by the CT, low-T DR, and RR rate coefficients. 

Taken together, the atomic data uncertainties alone can result in errors of nearly a factor of 
two (0.3 dex) in the derived Se abundances. In the case of Kr, if only [Kr III] is observed (as 
for most objects in our full sample), the abundance errors can be as large as 0.26 dex from 
the atomic data uncertainties. However, since we have empirically adjusted the Kr atomic 
data, the abundances are not likely to be so sensitive to atomic data uncertainties. 

This exercise reveals the sensitivity of our analysis to the unknown Se and Kr atomic 
data. Without a dedicated effort to calculate and experimentally measure these atomic data, 
uncertainties in the ionization balance for Se and Kr will not allow their abundances to be 
derived to accuracies of better than 0.2-0.3 dex. In fact, when other sources of error (e.g., 
observational uncertainties) are taken into account, their abundances may be considerably 
more uncertain. 



3.2. XSTAR Models 

In order to test the efficacy of the Se and Kr ionization balance and abundance results 
of our Cloudy models, we have modeled the ten objects listed in Table 1 with the pho- 
toionization code XSTAR (Kallman & Bautista 2001; Bautista & Kallman 2001). This code 
provides an independent test of our Cloudy results. 

XSTAR is a ID photoionization code, like Cloudy, but does not incorporate dust physics. 
Therefore, we used the same model assumptions (spherical geometry, constant density, no 
dust) as for the Cloudy models. Since XSTAR models are significantly more computationally 
demanding than Cloudy, we did not attempt to fit the observed spectra by optimizing the 
input parameters. Instead, we used the best-fit Cloudy values as input parameters to these 
models, although it occasionally was necessary to adjust the outer radii slightly so that the 
models were not too optically thick or thin. 

We find that the XSTAR models reproduce the predicted Cloudy spectrum for each PN 
reasonably well in all cases, with the intensities agreeing to within 30% for most lines. In 
particular, the [Se IV] 2.287 fim line intensities agree with the Cloudy results to within 30%, 
and often better, for all of the modeled PNc. The predicted Kr line intensities show larger 
differences from the Cloudy results, sometimes by as much as 50%. However, no systematic 
trends were apparent for the models with discrepant Kr line strengths. We therefore conclude 
that the Se and Kr abundance errors introduced by using different photoionization codes 
(<30-50%) are less than those of other sources of uncertainty (line intensities. Teg, densities, 
and atomic data errors). Our ionization balance and abundance results for Se and Kr thus 



-14- 



are not strongly model-dependent. 

4. IONIZATION CORRECTION FACTORS FOR SE AND KR 

The primary goal of this study is to derive widely applicable recipes to correct for the 
abundances of unobserved ionization stages of Se and Kr, so that their elemental abundances 
may be determined. For the vast majority of the PNe in our full sample (120 objects), we 
have observed only one ionization stage of each of these elements: Se^+ and Kr^"^. 

Therefore, we have constructed a grid of Cloudy models spanning a large range of phys- 
ical parameter space, to search for correlations between the fractional abundances of Se^+ 
and Kr++ and ionic fractions of more commonly detected elements. In particular, the grid in- 
cludes models for various values of Tea, nu, and ionization parameter U — Q(li)/{4:7rRf^nB.c), 
where (5(H) is the hydrogen- ionizing photon density, i?in the inner radius of the nebula, and 
c the speed of light. Given the density, the ionization parameter essentially determines the 
ionizing luminosity of the central star. 

We list in Table 10 the full set of physical parameters which were varied in our grid of 
models. A model was run for every combination of T^s, rin, and log([/) in Table 10, for a 
total of 3761 models. 

In our grid of models we have ignored dust, and assumed the default PN abundances 
in Cloudy (from Aller & Czyzak 1983 and Khromov 1989). We have tested the sensitivity 
of our results to metallicity and the presence of dust, which we discuss in the subsections 
below. All models have spherical geometry, constant density, an inner radius of 10~^'^ pc, 
and extend to an outer radius where Tc falls below 4000 K. Ranch (2003) NLTE stellar 
atmospheres with log{g) = 6.0 and solar abundances were used in all models. 

For each model, we extracted the fractional abundances (averaged over volume) of all 
ions of He, C, N, O, and Ne, and the first 10 ions of S, CI, Ar, Sc, and Kr, and compared 
the Se and Kr fractional abundances with several ionic fractions (or combinations thereof) 
of these other elements. In Figure 1, we plot the best correlations between the observed Sc 
and Kr ions and other ionic ratios. We also include correlations with (Kr++-|- Kr^+)/Kr, for 
completeness. In the six panels of Figure 1, the fractional abundances from each model in 
our grid are represented as dots. In some of the panels, e.g., Se^+/Se vs. a family of 

curves is seen among the plotted ionic fractions. Each curve corresponds to a single value of 
Teff, indicating the correlation's sensitivity to the nebular gas temperature (which generally 
increases with Teflf)- 
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We have performed fits to the ionic fraction correlations, using a minimization routine 
written in IDL, and considering polynomial and exponential functional forms. These are 
displayed in Figure 1 as solid lines. The inverses of these fitting functions correspond to 
ICF formulae which can be used to convert Se and Kr ionic abundances into elemental 
abundances. 

In particular, when [Kr III] is the only Kr ion observed, the Kr ICF may be derived 
using the fractional abundance of either S"'""'" or Ar"*""*", which are commonly detected in 
optical spectra of PNe. The ICFs are: 

ICF(Kr) = Kr/Kr++ = (-0.009205 + 0.3098x + 0.0007978e^-2'^^^)-\ (1) 

X = Ar++/Ar > 0.027; 

and 

ICF(Kr) = Kr/Kr++ = (-0.3817 + 0.3796e^-°^22/)-i^ (2) 
y = S++/S > 0.0051. 

The lower hmits to x and y denote the range of validity for these ICFs (below these limits, the 
ICFs become negative). Of these two ICFs, we recommend Equation (1) over Equation (2) 
since more ions of Ar (Ar"*""*", Ar^"*", and Ar^"*") are typically observed in optical spectra of 
PNc than of S and S"'"'^). This leads to a more accurate Ar abundance from optical 
data than that of S. In fact, Henry ct al. (2004) have found that S abundances in PNe may 
be inaccurate when determined from optical data and model-derived ICFs, and explain this 
behavior by the inability of photoionization models to accurately compute the abundance of 
S^"*", which can only be observed with IR spectroscopy. 

Nevertheless, both correlations (top panels of Figure 1) are quite strong, and if the 
Ar and S abundances are well-determined, both should be reliable if our models accurately 
reproduce the Kr ionization equilibrium. In Paper 11, we make use of both of these equations. 

On the other hand, correlations between Se^"''/Se and ionic fractions of other elements 
are not as well-defined. This will lead to less accurate ICFs, and hence less accurate Se 
abundance determinations. In the middle left panel of Figure 1, we display the fit to the 
best correlation we found, between Se^+/Se and 0"'""''/0. This corresponds to 

ICF(Se) = Se/Se=^+ = (-0.1572 - 0.3532z^^-^^ + 0.153e^-^^^^)-\ (3) 

z = 0++/0 > 0.01626. 

We have also derived ICFs to be used when both [Kr III] and [Kr IV] are observed. 
Although we do not utihze these ICFs in our study, we present them for completeness. The 
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ICF that is least sensitive to variations in rie, T^s, and U is 

ICF(Kr) = Kr/(Kr++ + Kr^+j = (-0.04996 + 0.7041a; + QM79x^)-\ (4) 

X = (C1++ + Cl^+)/Cl > 0.0686. 

However, in order to detect [CI IV] , it is necessary to have obtained a spectrum extending 
to at least 7550 A. For data which do not extend to such long wavelengths, we offer the 
following ICFs: 

ICF(Kr) = Kr/(Kr++ + Kr^+j = (-0.1549 + 0.6458y + 0.0048e^-^^^^)-\ (5) 

y = (Ar++ + Ar^+)/Ar > 0.2193; 
ICF(Kr) = Kr/(Kr++ + Kr^+) = (-0.2074 + 0.9834z + 0.000513e^-^^2^)-\ (6) 

(0+ + 0++)/0 > 0.21. 

Both of these correlations exhibit large dispersion and a non-trivial dependence on T^s and 
log(t/), and are less accurate than Equation (4) unless y and z are large (the dispersion 
about the fit becomes small in this case; see the lower panels of Figure 1). 

4.1. Dependence of ICFs on Metallicity 

To derive Equations (l)-(6), we have used the default PN abundances of Cloudy. How- 
ever, some of the PNe of our sample may have substantially different abundances. In order 
to test the robustness of our ICFs to the assumed metallicity, we have run the same grid of 
models for low metallicity PNe by dividing the default Cloudy PN abundances by 10 for all 
elements heavier than He, and using NLTE stellar atmospheres with halo abundances (Ranch 
2003). Note that these abundances were not chosen to represent specific PNe, but are simply 
used to illustrate that our derived ICFs are not strongly dependent on the adopted nebular 
abundances. 

In Figure 2, wc show the Sc and Kr ionic fraction correlations discussed above when 
low metallicity abundances are adopted. We overplot the fits to the correlations of Figure 1 
as solid lines. As can be seen from Figure 2, adopting lower metal abundances does not 
significantly affect the correlations we found in the previous section. In particular, the 
correlations described by Equations (l)-(4) for high metallicities are remarkably similar to 
those at low metaUicities (although slightly less Se^"*" is produced at high fractions in 
the low metallicity models). The only significant discrepancies in the correlations are for 
those in the lower panels of Figures 1 and 2, where we find higher (Kr"*""*"-!- Kr^"*") fractional 
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abundances in our low metallicity models for low (Ar+''"+ Ar^+)/Ar and (0++ 0++)/0. 
This signifies a deterioration in the accuracy of Equations (5) and (6) for low metallicity, 
low excitation PNe. However, these correlations exhibit a large dispersion even at higher 
metallicities, and are not as reliable as the correlation of Equation (4), which applies to low 
metallicity PNe. We conclude that adopting a lower metallicity does not significantly affect 
the ICFs of Equations (l)-(4). 

4.2. Dependence of ICFs on Dust 

To test the dependence of our ICFs on the inclusion of dust physics, we have run our 
grid of models with C-rich grains (the "grain Orion graphite" command in Cloudy) and 
0-rich dust (the "grain Orion silicate" command), with grain properties from Baldwin et al. 
(1991) and van Hoof et al. (2004). 

In Figure 3, we display the ionic fraction correlations for Cloudy models including C-rich 
dust. While some minor differences in the fractional abundances of these ions arc apparent 
when compared to Figure 1, the ICFs derived from models without dust (overplotted as solid 
lines) also apply to models that include dust. 

In our 0-rich dust models, the silicate grains began sublimating for high values of log(f/) 
(> 1.0), and hence we were not able to run our full grid. However, the correlation plots (not 
shown) are similar to those of the C-rich dust models, and again no significant differences 
are seen compared to the dust-free models. We therefore conclude that the presence of dust 
does not affect our derived ICFs. 



4.3. Compcirison of ICF- and Model-Derived Se and Kr Abundances 

In this section, we compare the Se and Kr abundances determined from our ICFs to 
those found from our models of the ten PNe listed in Table 1 (§3.1). Since our interest is 
primarily in determining Se and Kr abundances from their NIR lines, we only consider the 
ICFs given by Equations (l)-(3). 

For each of the ten PNe we have modeled, we have taken ionic and elemental abundances 
needed for our ICFs (Ar"'""''/Ar, S^"''/S, and 0+"'"/0) from the literature (see Table 1). To 
properly account for uncertainties in the derived ICFs, we consider errors in these ionic and 
elemental abundances. Unfortunately, such an error analysis has not been performed in any 
of the abundance references we use. Therefore, we adopt 20% uncertainties in the O, S, and 
Ar ionic and elemental abundances (see Paper II). We propagate these errors into our ICFs, 
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as well as the standard deviations of our fits to the ionic fraction correlations in Figure 1. 

To determine the Kr++ and Se^+ ionic abundances, we have utilized 5-level and 2-level 
model atoms, respectively, with transition probabilities and effective collision strengths from 
the references listed in Table Bl. In these calculations, we have used the T'e([0 III]) and Ue 
derived in the optical data references. Uncertainties of ±1000 K in Te and 20% in Ue were 
assigned when not reported in these references, and propagated into the Se^"^ and Kr++ 
abundance errors. 

We present the results of our calculations for Kr in Table 11 and Se in Table 12. In 
these tables, the ICFs derived from Equations (l)-(3) are given, followed by the abundances 
determined using our ICFs and photoionization models. Despite the considerable uncer- 
tainties we estimate for the ICFs, we find gratifying agreement between the ICF and model 
abundance determinations. In all cases, the abundances agree within the errors, and with a 
couple of exceptions agree to within 0.15 dex (~40%) for Kr and 0.1 dex (~25%) for Se. This 
also illustrates that the simplifying assumptions we made in our models do not significantly 
affect the derived Se and Kr abundances. 

The largest discrepancies in the Kr abundances occur for objects in which only 
[Kr III] A6826.7 was detected, while 2.199 /xm was unobserved. As discussed in Appendix A. 3, 
A6826.7 can be blended with He I A6827.9 and telluric OH emission in low resolution spectra. 
Therefore, the intensity of [Kr III] A6826.7 may be quite uncertain if these contaminants are 
not properly removed. 

Another major source of error in our ICF-derived abundances stems from the ICF itself. 
In some cases, the uncertainties in the ICF approach a factor of two. Furthermore, for small 
or Ar++ fractions, the ICF becomes large, and may conceivably be more uncertain 
than we have estimated. We befieve that the uncertainties in the ICFs and [Kr III] A6826.7 
intensities are able to explain the largest discrepancies between our ICF- and model-derived 
Se and Kr abundances, namely for NGC 6741, NGC 6826, and NGC 7662. Nevertheless, the 
overall agreement between the Se and Kr abundances derived with these different methods 
is encouraging, and validates the ICFs we have derived. 

5. CONCLUDING REMARKS 

We have presented results of a photoionization model study designed to derive correction 
factors for the abundances of unobserved Se and Kr ions. These corrections are needed to 
determine elemental abundances from our survey of [Kr III] 2.199 and [Se IV] 2.287 fira in 
120 Galactic PNe, the first large-scale survey of n-capture elements in PNe. We have added 
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Se and Kr to the atomic databases of the photoionization codes Cloudy and XSTAR. While 
the atomic data concerning energy transitions for the observed Se and Kr ions are known 
(i.e., transition probabilities and effective collision strengths), much of the data governing the 
ionization equilibrium of these elements have not been determined. We have approximated 
the rate coefficients for these processes as outlined in Appendix B. 

When deriving ICFs, it is critical to accurately solve for the ionization balance of the 
element in question. To this end, we calculated detailed models of ten PNe exhibiting emis- 
sion lines from multiple ionization states of Kr. The Kr atomic data were then empirically 
adjusted until the emission line intensities from the various Kr ions were adequately repro- 
duced by our models. Unfortunately, no other Se ions have been clearly detected in PNe, 
and thus we could not empirically correct the uncertain Se atomic data we adopted. This 
makes the Se abundances we derive from our survey inherently more uncertain than those 
of Kr. 

Using the adjusted Kr atomic data, we computed a grid of Cloudy models spanning a 
large range of T^s, riu, and ionization parameter U, which encompass the values of these 
parameters hkely to be encountered in most PNe. Prom these grids, we extracted the frac- 
tional abundances of the observed Se and Kr ions, and compared them to the ionic fractions 
of commonly detected elements. We computed fits to the best correlations we found, which 
serve as ICFs that can be used to determine Se and Kr elemental abundances when only one 
or two ions of these species are observed. The ICF formulae are given in Equations (l)-(6). 

We have compared the Se and Kr abundances derived with the ICF method with those 
determined from our photoionization models of ten PNe. In most cases, the abundances 
agree to within 25-40%. The few discrepant abundances serve to illustrate the magnitude 
of errors in the ICFs. When uncertainties in the ionic and elemental abundances used in the 
ICFs arc taken into account, the Se and Kr ICFs are uncertain by at least 20% and 40%, 
respectively, and the errors may approach a factor of two in some cases. It should be noted 
that the ten PNe we have considered are among the brightest and most well-studied PNe in 
the Northern Hemisphere. For fainter objects in our NIR survey, the ICFs may be even more 
uncertain. Therefore, when a realistic and rigorous error analysis is applied, our derived Se 
and Kr abundances are unlikely to be determined with an accuracy of better than 50%. 

A more serious factor in the accuracy of our Se and Kr abundance determinations is 
the uncertain atomic data governing the ionization equilibrium of these elements (PI cross- 
sections and rate coefficients for RR, DR, and CT). We have run Monte Carlo simulations 
to test the importance of these uncertainties in our models of four PNe, allowing the rate 
coefficients of these atomic processes to vary within a Gaussian distribution of FWHM 
0.5 dex around the adopted values, or a factor of three (1-cr) uncertainty. In particular. 
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uncertainties in the PI cross-sections are the most significant, and can lead to errors of up 
to 0.2-0.25 dex (~ 60-80%) in the predicted fine intensities of [Se IV], [Kr III], [Kr IV], 
and [Kr V]. Uncertainties in the CT rate coefficients are the next most important, especially 
for the lower ionization states, and can result in errors of up to 0.15 dex (40%) in the 
[Kr III] intensities. The rate coefficients for low-T DR are also significant, and can lead to 
errors of 0.1-0.2 dex (25-60%) in the predicted intensities of high ionization Se and Kr lines, 
particularly for low excitation objects. The uncertainties in the predicted line intensities 
correspond to uncertainties in the abundances of the parent ions, which propagate into the 
elemental abundance errors. Taken together, the atomic data uncertainties alone can result 
in Se abundance errors of 0.3 dex (a factor of two), and Kr abundance uncertainties of up 
to 0.2-0.25 dex. 

Therefore, calculations and measurements of PI cross-sections and RR, DR, and CT rate 
coefficients for Se and Kr ions are needed to improve the accuracy to which the abundances of 
these elements may be determined in PNc. Furthermore, most other n-capture elements also 
lack such atomic data, and this is a serious impediment to studying s-process enrichments in 
PNc. Other n-capture elements, including Br {Z = 35), Rb {Z = 37), and Xe {Z = 54) have 
been detected in NGC 7027 (PB94; Zhang et al. 2005) and other bright PNe (Sharpee et al. 
2003; Liu et al. 2004a; Zhang et al. 2006). Accurately determining the abundances of these 
elements, particularly Xe, in conjunction with Se and/or Kr would provide crucial details 
of the physical conditions under which s-process nucleosynthesis operates in PN progenitor 
stars. 

An investigation of the unknown Se and Kr atomic data is difficult and time-consuming, 
and was not attempted for our current study. However, one of us (NCS) has begun a 
laboratory astrophysics study of Se, Kr, and Xe, which will provide much of the unknown 
atomic data for these elements. This study will utilize both theoretical calculations and 
experimental measurements to determine the data. In fact, the PI cross-sections of Se"*" and 
Se"'"^ near threshold have already been measured at the Advanced Light Source in Berkeley, 
CA, and are discussed in Sterling (2006b). 

While the atomic data uncertainties must be acknowledged as an important source of 
error in our Se and Kr abundance determinations at this stage, our model-derived ICFs are 
as accurate as is currently possible. In Paper II, we will apply these ICFs to our full sample 
of PNe to derive Se and Kr abundances or upper limits. The large fraction of PNc in which 
the Se and/or Kr NIR lines have been detected (81 of 120 objects) enables us to determine n- 
capture element abundances in a large sample of Galactic PNe for the first time; indeed, our 
survey will increase the number of PNe with known n-capture element abundances by more 
than a factor of ten. With this data, we provide in Paper II the first broad characterization 
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of s-process enrichments in PNe, which reveals their dredge-up histories and importance to 
the evolution of trans-iron elements in the Galaxy. 
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A. APPENDIX: SE AND KR LINE IDENTIFICATIONS 

In this appendix, we discuss the identifications of the optical and NIR Kr and Se emission 
lines listed in Table 2. These lines have been associated with several different species in the 
literature, if identified at all. For example, [Kr IV] A5346.0 has previously been identified 
as G III (Liu et al. 2004a) and S II (Peimbert et al. 2004; Hyung et al. 1994a, 1995), while 
[Kr IV] A5867.7 has been associated with a permitted Al II triplet (Sharpee et al. 2003) as 
well as He II Pfund 29. 

Gonsidering the confusion in the literature regarding the identity of these lines, we 
deemed it worthwhile to investigate all possible identifications for the putative Se and Kr 
fines. We have determined that, at least for the 10 PNe we have modeled in detail, the 
evidence wc have accumulated strongly indicates that the lines listed in Table 2 are indeed 
from Se and Kr. 

In the following, we consider each line separately. To illustrate our identification method, 
we present our investigation of A5346.0 in detail. For the other features, we provide less detail, 
but discuss important line blends and how the contaminants were removed. 
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A.l. The A5346.0 Feature 

We have used the Atomic Line List v2.0^ compiled by P. A. M. van Hoof to search for 
ahernate identifications of [Kr IV] A5346.0. We searched a region ±0.5 A around the central 
wavelength, and considered all elements through Z = 30. We excluded forbidden transitions 
between levels higher than 2 eV above the ground state, as well as intercombination lines of 
elements other than C, N, and O. Furthermore, we did not consider permitted transitions of 
rare or highly refractory elements (e.g., Cr, Cu, Mg, Mn). The remaining possible identifi- 
cations are listed in Table Al, along with the energies of the upper and lower levels of the 
transition, and the wavelengths of other lines from the same multiplet. This list includes 
the tentative S II and C III identifications that have appeared in the literature (Hyung et al. 
1994a, 1995; Liu et al. 2004a; Peimbert et al. 2004). 

With the above constraints, all of the possible alternative identifications of A5346.0 are 
permitted lines. If any of these identifications are correct, other lines from the same multiplet 
should be seen. Of course, some of these may be undetected if they are relatively weak, and 
for many of the species in Table Al, the relative strengths of other multiplet members are 
unknown. Nevertheless, it is noteworthy that we do not find any of the other multiplet 
members in the optical spectra of the ten modeled PNe. 

In particular, we can dismiss the identification as C III A5345.9. According to Wiese et 
al. (1996, hereafter WFD96), in LTE C III A5359.9 should be nearly twice the strength of 
A5345.9. Although this region of the spectra is free from other features, this line is undetected 
in all of the optical spectra we have utilized. C III A5337.4, half the strength of A5345.9 in 
LTE (WFD96), is also unseen. These lines are also undetected in echelle spectra of eight 
of the modeled PNe (all but NGC 6826 and NGC 6886) obtained by Dinerstein et al. (in 
preparation; see also French et al. 2000), despite the presence of A5346.0 in each. 

Analogously, N I A5346.5 is the weakest member of its triplet in LTE (WFD96), but the 
strongest line, A5310.4, is undetected in the optical spectra of all 10 PNe. 

While the relative strengths of multiplet members for other lines in Table Al are not 
known (to our knowledge), we consider their non-detection to be evidence that these species 
do not contribute significantly to the fiux of A5346.0. The only plausible identification that 
remains is [Kr IV] A5346.0, which is strengthened by the presence of another [Kr IV] line at 
5867.7 A. 



■^P. A. M. van Hoof (1999) , |http: / / www.pa.uky.edu/ ~peter /atomic / 
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A.2. The A5867.7 Feature 

As with [Kr IV] A5346.0, we compiled all possible alternate identifications of 
[Kr IV] A5867.7 with the Atomic Line List. We excluded Al II AA5867.6, 5867.8, 5867.9 
(Sharpee et al. 2003) based on the non-detection of their multiplet members at AA5853.8 
and 5861.7 (although the latter may be lost in a blend with [Mn V] A5861.0 in some ob- 
jects). Other possible contaminants were ruled out using the same argument, with one 
exception. 

He II A5869.0 (Pfund 29) is a potential contaminant to [Kr IV] 5867.7, especially in low 
resolution spectra (such as that of Liu et al. 2004a,b and Zhang et al. 2005). Therefore, we 
have subtracted the contribution of this hne to the reported A5867.7 intensities in the optical 
data references, by using the intensities of nearby isolated He II Pfund series members and the 
Fortran program of Storey & Hummer (1995). Note that we perform this subtraction even in 
spectra with sufficiently high resolution to resolve these features if the He II intensity was not 
given separately from that of [Kr IV]. In these calculations, we used the [O III] temperatures 
and averaged electron densities derived from the optical spectra. The remaining fiux, which 
we attribute to [Kr IV], is given in Table 2. This correction reduces /(A5867.7) by as little 
as 8% (NGC 6790) to up to 40% (NGC 6886), and hence [Kr IV] is the primary contributor 
to this feature in all the PNe we have modeled. 



A.3. The A6826.7 Feature 

Baluteau et al. (1995) listed two potential contaminants to [Kr III] A6826.7: He I A6827.9 
and C I A6828.1. These transitions are the most likely alternate identifications for this line, 
and other possibilities have been excluded as in §A.l. 

C I A6828.1 is a singlet; however, we deem this an unlikely contaminant since other lines 
from the same upper level (AA6294.1, 6656.5, 7366.8) are not seen in the optical spectra of 
the 10 PNe we have modeled. 

On the other hand. He I A6827.9 is likely to contribute to this line in low resolution 
spectra, and other lines from the S^S^^P series (Pequignot & Baluteau 1988; Smits 1991) 
have been detected in the modeled PNe. Using the intensities of neighboring lines in this 
series and the Smits (1991) tables of theoretical He I hne intensities at Te — 10^ K and 
He — 10^ cm~^, we have removed the contribution of He I A6827.9 to the A6826.7 feature. In 
some cases (NGC 6790, NGC 6826, NGC 6884, NGC 7662) this correction can be significant 
(~40%), while for the other PNe of our sample, it reduces /(A6826.7) by less than 20%. We 
attribute the remaining flux (fisted in Table 2) to [Kr III] A6826.7. 
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In addition to nebular contaminants, [Kr III] A6826.7 can also be affected by telluric OH 
emission at A6827.5 (Hanuschik 2003; Osterbrock et al. 1996). For this reason, the intensity 
of [Kr III] A6826.7 from low resolution spectra (e.g., Liu et al. 2004a,b; Zhang et al. 2005) 
is likely to be more uncertain than that of 2.199 fira. This telluric feature may also affect 
the A6826.7 flux in two other PNe exhibiting [Kr IV] emission, NGC 2022 and NGC 6818 
(Tsamis et al. 2003). These spectra have resolutions of 4.5 A FWHM at this wavelength, 
which is sufficiently low that blending with telluric features may be important. When we 
attempted to model the Kr features in these two PNe, we found that the [Kr III] A6826.7 
intensity was far too strong (by more than an order of magnitude) to be consistently fit with 
the [Kr IV] lines. Since no He I 3^S-n^P series lines were detected in these two PNe, the 
discrepancy may be due to incomplete removal of the telluric emission near A6826.7, or an 
instrumental artifact. We did not attempt to model these two PNe further. 



A.4. The 2.199 /^m Feature 

Identified as [Kr III] by Dinerstein (2001), this line had long been a mystery to infrared 
astronomers. It was first detected in a PN by Treffers et al. (1976), and Geballe et al. 
(1991) unsuccessfully attempted to identify it (and the 2.287 /xm feature) with transitions of 
elements from H to Ca and Fe. While H2 3-2 S{3) 2.201 fim possibly contributes to the flux 
of this line in some PNe, the detection of the 2.199 fim feature in PNe that do not exhibit 
molecular hydrogen indicates that it has an atomic nature, as noted by Geballe et al. 

Cox et al. (2002) identified this line as N I 2.199 ^m in NGC 7027, based on the detection 
of other N I lines at 2.194 and 2.148 fim. However, we do not detect the latter N I lines in 
any of the 103 PNe we observed, nor are they evident in any of the literature spectra we 
have utilized, including those of NGC 7027 (Geballe et al. 1991; Hora et al. 1999; Lumsden 
et al. 2001). Examining Figure 1 of Cox et al. (2002), one sees that N I 2.194 /xm is much 
weaker than the 2.199 /xm fine, and N I 2.148 /xm lies in a region of their spectrum which 
is extremely noisy, and hence its flux is quite uncertain. This suggests that N I is a minor 
contaminant to [Kr III] 2.199 fxm even in NGC 7027; the absence of NIR N I features in all 
of the other PNe in our sample indicates that the N I contribution to this line is negligible 
in most cases, and we do not consider it further. 

Similarly, we flnd that other transitions of H through Zn within 0.005 /im of the 
2.199 fim line are unlikely contaminants, based on the non-detections of permitted line 
multiplet members, and the lack of forbidden transitions between low-lying energy levels in 
this wavelength region. Therefore, the only important contaminant to [Kr III] 2.199 //m is 
H2 3-2 S{3) 2.201 urn. 
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Of the ten objects considered in this paper, four exhibit NIR H2 emission (IC 5117, 
NGC 6741, NGC 6886, and NGC 7027). The contribution of H2 to the 2.199 fim flux 
depends on the H2 excitation mechanism. For example, models of Black & van Dishoeck 
(1987) show that the strengths of H2 3-2 lines in the K band are negligible if the H2 is 
thermally excited, but may be a significant fraction (30-40%) of the H2 1-0 >S'(0) 2.224 fim flux 
under fluorescent excitation. In principle, the H2 excitation mechanism can be ascertained 
from the flux ratio of H2 2-1 S{1) 2.248 //m relative to H2 1-0 S{0) 2.224 /xm, which we 
denote as F(2. 248/2. 224). The canonical fluorescent-excitation model. Model 14 of Black 
& van Dishoeck, predicts that F(2. 248/2. 224) = 1.22, while their thermal excitation model 
(S2) predicts F(2.248/2.224) = 0.38. The four PNe considered in this paper that exhibit 
H2 emission all have F(2. 248/2. 224) < 0.65, indicating that the observed H2 is collisionally 
excited, and hence contamination of the [Kr III] flux (as well as [Se IV]) by H2 is not expected 
to be important in these objects. 



A.5. The 2.287 //m Feature 

First identified as [Se IV] by Dinerstein (2001), this hue also may be contaminated by 
molecular hydrogen (H2 3-2 S{2) 2.287 //m), but in §A.4 we show this to be unimportant 
in the ten objects we consider in this paper. Our search for additional contaminants of this 
line revealed no other likely identifications, in agreement with the conclusions of Geballe et 
al. (1991). We searched for transitions from H through Zn within 0.005 /^.m of 2.287 /xm, 
and the only likely candidates were permitted lines for which other members of the multiplet 
were in all cases undetected. 



A.6. The AA6256.1, 8243.4 Features 

These two wavelengths correspond to the strongest optical [Kr V] lines. Since A6256.1 
was detected in only one PN in our sample (NGG 7027) and A8243.4 not at all, we have not at- 
tempted detailed investigations of the identities of these features. We only note that one must 
use caution in identifying A6256.1 with [Kr V] due to the nearby fines C II AA6256.2, 6256.5, 
and that [Kr V] 8243.4 is prone to blending with telluric emission, as well as N I A8242.4 
and O III A8244.1 (Baluteau et al. 1995). 
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B. APPENDIX: SE AND KR ATOMIC DATA 

In this appendix, we describe the adopted Se and Kr atomic data used in our Cloudy and 
XSTAR models. The atomic data governing the ionization equilibrium of these elements, 
and indeed all n-capture elements, are poorly known. In most cases, the photoionization (PI) 
cross-sections and rate coefficients for radiative recombination (RR), dielectronic recombina- 
tion (DR), and charge transfer (CT) have not been calculated or experimentally measured 
for Se and Kr ions. For this reason, we have been forced to use approximate atomic data 
for these processes in order to model the ionization balance of Se and Kr. In the case of 
Kr, multiple ionization stages have been observed in the 10 PNe listed in Table 1, which 
allows us to adjust the uncertain Kr atomic data so that the line intensities are satisfactorily 
reproduced by our models. However, no other Se ions have been detected in PNe, and we 
are unable to empirically correct the Se atomic data. 

B.l. Energy Levels, Transition Probabilities, and Effective Collision Strengths 

Unlike the atomic data controlling the ionization balance of Se and Kr, emission line 
data has been determined for the observed ions of these elements. Transition probabilities 
Aij and effective collision strengths have been determined for Se IV, Kr III, Kr IV, and 
Kr V. In Table Bl we list the references from which these data are taken. 

All energy level data are taken from the NIST Atomic Spectra Databas^, as are the 
ionization potentials (IPs) with the exception of Kr^^^-Kr^^"^, whose IPs are taken from the 
calculations of Biemont et al. (1999). 

B.2. Photoionization Cross-Sections and Radiative Recombination Rate 

Coefficients 

PI cross-sections are unknown for all Se and Kr ions except Kr° (Richter et al. 2003) 
and Se° (Chen & Robicheaux 1994, who calculated the cross-section for only a small energy 
range near threshold). 

The classical Kramers (1923) formula, 

ath = 2^3-^-^a7ra^(£;Ryd/^th)(iVe/n), (Bl) 



''National Institute of Standards and Technology Atomic Spectra Database v3.0; see 
http://physics.nisLgov/PhysRefData/ASD/index.htnil 
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can be used to estimate the PI cross-section at the threshold energy (i.e., ionization po- 
tential). In this equation, ath is the threshold cross-section, a the fine-structure constant, 
Oq the Bohr radius, i?Ryd the Rydberg energy, i?th the threshold energy, iVg the number 
of valence electrons, and n the principal quantum number of the valence electrons. Gould 
(1978) has shown that this equation provides a reasonably good approximation (to within 
a factor of two) of the threshold PI cross-sections of first and second row elements. How- 
ever, Equation (Bl) does not provide information regarding the functional dependence of the 
cross-sections on photon energy. Therefore, we simply assumed that the PI cross-sections 
have power law forms above the threshold energies (Osterbrock 1989, pp. 34-36): 

a^E) = ath(^/^th)-'. (B2) 

Although (Jpi{E) is known for Kr° and (partially) for Se", for consistency we have used 
cross-sections from Equations (Bl) and (B2) for these species in Cloudy and XSTAR. With 
XSTAR, we have tested the sensitivity of our results to the adopted form of the Kr° and 
Se° PI cross-sections. We smoothed the Se° and Kr° cross-sections (Chen & Robicheaux 
1994; Richter et al. 2003) using the resonance-averaged photoionization (RAP) cross-section 
method of Bautista et al. (1998). For energies not considered by Chen & Robicheaux, we 
assumed (7pi(-E') is a power-law with index —2, as in Equation (B2). The Kr° and Se" PI 
cross-sections were entered into XSTAR as energy-cTp; pairs which are linearly interpolated 
to reproduce the cross-sections. As expected, we find that the predicted Se and Kr line 
intensities depend negligibly (< 2%) on the adopted Se° and Kr° PI cross-sections. The 
reason is that neutral Se and Kr are trace species in our models, and hence their fractional 
abundances do not significantly affect the derived elemental abundances. 

We have determined o"pi(i?) for all terms in the ground principal quantum number shell 
of each Se and Kr ion with Equations (Bl) and (B2), and used the Milne relation (e.g., 
Osterbrock & Ferland 2006, p. 401) to determine the RR rate coeflicient A^,- into this shell. 
For higher principal quantum number shells, we adopted hydrogenic recombination rate 
coefficients (Seaton 1959). Adding the A^^ for the ground and higher principal quantum 
number shells of an ion gives the total ^4^. The A„ of high ions of Se and Kr (greater than 
5-6 times ionized) are underestimated, due to incomplete energy level data for these species 
(the lack of energy level data in the ground principal quantum number shell prevents the 
calculation of RR rate coefficients into these states). However, since these ionization states 
are usually sparsely populated in photoionized nebulae, this should have a negligible effect 
on our results. 

The dependence of the RR rate coefficient on temperature is often approximated as a 
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power law (e.g., Aldrovandi & Pequignot 1973), 

ar{T) = A„r^ (B3) 

where t is the electron temperature in units of 10^ K. We adopt rj = 0.7, which has been 
found to be a good approximation even for heavy element ions (e.g., Mazzitelli & Mattioli 
2002). 

The Se and Kr PI cross-sections and RR rate coefficients we use in Cloudy and XSTAR 
are given by Equations (B1)-(B3), with the fit parameters listed in Table B2. 

B.3. Charge Transfer Rate CoefRcients 

CT rate coefficients are not known for any ions of Se and Kr. The heaviest elements 
with reliable CT data are Ni and Cu (Kingdon & Ferland 1996)|f|. Therefore, we assume that 
the rate coefficients of the first four ions of Se and Kr are the same as the averaged rate 
coefficients of similar charge states of Ni and Cu. 

For higher Se and Kr ions, we adopt the formula of Ferland et al. (1997): 

act = 1-92 X 10-^C (cm^ s"^), (B4) 
where is the rate coefficient, and ( is the ion charge. 

B.4. Dielectronic Recombination Rate CoefRcients 

DR is the dominant recombination process for many species in photoionized nebulae. 
There are essentially two different types of DR which can be important in ionized plas- 
mas: high-temperature DR, where core relaxation stabilizes the ion, and low-temperature 
DR, where the recombining electron is captured into a low-lying resonance state (Dopita & 
Sutherland 2003, p. 108). 

In photoionized nebulae, low-T DR dominates high-T DR in almost all cases. However, 
low-T DR rate coefficients are unknown for most third row and higher elements, because the 



^While these authors derived CT data for Zn, the rate coefficients of Zn'^^ and Zn'*+ are iilceiy to be 
underestimates due to the incomplete energy level data of these ions (since electrons are expected to be 
captured into energy levels higher than the limit of the published levels). 



-29- 



energies of the low-lying autoionizing resonance states have not been experimentally mea- 
sured (Ferland et al. 1998; Ferland 2003). The unknown energies of these states precludes 
theoretical calculations of the low-T DR rate coefficients for heavy elements. In principle, 
the rate coefficients can be determined experimentally with electron-ion beam merging tech- 
niques such as those of Savin et al. (1999, 2006), but it is difficult to use these methods 
to measure DR rate coefficients for the low charge-to-mass species found in photoionized 
nebulae, and to date no such measurements have been performed. 

In Cloudy, low-T DR is approximated for all third row and higher elements in a manner 
depending on the ionization stage. For the first four ions, the rate coefficients are assumed 
to be equal to the averaged rate coefficients of the same charge states of C, N, and O. For 
higher ions, the rate coefficient is given by 

a^ow-T ^ 3 Q ^ iqO.ic-12 ^^^3 

where ( is the ion charge. We have adopted these approximations for Se and Kr in Cloudy, 
and added them to the atomic database of XSTAR for all third row and heavier elements. 

Although high-T DR is less important in photoionized nebulae, rate coefficients for this 
process have been determined for all Kr ions by Fournier et al. (2000). We fit these data 
with the formula of Aldrovandi & Pequignot (1973, their equation 5): 

cyT-"" = ^diT-i-^e-^o/^^(l + i?die-^^/^0, (B6) 

where is the electron temperature (in K), and Adi, Tq, and Ti are the fitting parameters 
(listed in Table B3). The high-T DR rate coefficients for Se are unknown, and hence we 
assumed the same rate coefficients as Kr ions of the same charge. 



B.5. Collisional Ionization Cross-Sections 

Like high-T DR, collisional ionization (CI) plays a minor role in most photoionized 
nebulae, since the plasma is at too low a temperature (~1 eV) to ionize atoms. For com- 
pleteness, we have added the Kr CI data of Loch et al. (2002)0 to Cloudy and XSTAR. We 
fit the data with the formula of Voronov (1997): 

< ^c^^ >= ^ci%^^^f/^e-^ (B7) 



http: / / www-cfadc.phy.ornl.gov /data^nd_codes / aurost /aurost Joniz /kr-iso- nuclear .html 
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where < a^v > is the rate coefficient, U — Eth/T^, and A^, P, X, and K are tfie fit 
parameters (given in Table B3). Tlie CI rate coefficients are not known for Se, and lience 
we assume Se ions fiave tfie same rate coefficients as Kr ions of tfie same cfiarge. 

REFERENCES 

Aldrovandi, S. M. V., & Pequignot, D. 1973, A&A, 25, 137 

Aller, L. H., & Czyzak, S. J. 1983, ApJS, 51, 211 

Aller, L. H., Hyung, S., & Feibelman, W. A. 1996, PASP, 108, 488 

Asplund, M., Grevesse, N., & Sauval, A. J. 2005, in "Cosmic Abundances as Records of 
Stellar Evolution and Nucleosynthesis" , ASP Conf. Ser. 336, eds. T. G. Barnes III & 
F. N. Bash (San Francisco: ASP), 25 

Baldwin, J. A., Ferland, G. J., Martin, P. G., Corbin, M. R., Cota, S. A., Peterson, B. M., 
& Slettebak, A. 1991, ApJ, 374, 580 

Baluteau, J.-R, Zavagno, A., Morisset, C, & Pequignot, D. 1995, A&A, 303, 175 

Bautista, M. A., & Kallman, T. R. 2001, ApJS, 134, 139 

Bautista, M. A., Romano, P., & Pradhan, A. K. 1998, ApJS, 118, 259 

Biemont, E., Fremat, Y., & Quinet, P. 1999, Atom. Data Nuc. Data Tab., 71, 117 

Biemont, E., & Hansen, J. E. 1986a, Phys. Scr, 33, 117 

Biemont, E., & Hansen, J. E. 1986b, Phys. Scr, 34, 116 

Biemont, E., & Hansen, J. E. 1987, Nuc. Inst. Meth. Phys. Res. B, 23, 274 

Black, J. H., & van Dishoeck, E. F. 1987, ApJ, 322, 412 

Busso, M., Gallino, R., & Wasserburg, G. J. 1999, ARA&A, 37, 239 (BGW99) 

Busso, M., Gallino, R., Lambert, D. L., Travaglio, C, & Smith, V. V. 2001, ApJ, 557, 802 

Cahn, J. H., Kaler, J. B., & Stanghellini, L. 1992, A&AS, 94, 399 

Cardelh, J. A., Federman, S. R., Lambert, D. L., & Theodosiou, C. E. 1993, ApJ, 416, L41 
Chen, C.-T., & Robicheaux, F. 1994, Phys. Rev. A, 50, 3968 



-31 - 



Cox, P., Huggins, P. J., Maillard, J.-R, Habart, E., Morisset, C, Bachiller, R., & Forveille, 
T. 2002, A&A, 384, 603 

Dinerstein, H. L. 2001, ApJ, 550, L223 

Dopita, M. A., & Sutherland, R. S. 2000, ApJ, 539, 742 

Dopita, M. A., & Sutherland, R. S. 2003, Astrophysics of the Diffuse Universe (Berlin: 
Springer) 

Ferland, G. J. 2003, ARA&A, 41, 517 

Ferland, G. J., Korista, K. T., Verner, D. A., & Dalgarno, A. 1997, ApJ, 481, L115 

Ferland, G. J. Ferland, Korista, K. T., Verner, D. A., Ferguson, J. W., Kingdon, J. B., & 
Verner, E. M. 1998, PASP, 110, 761 

Fournier, K. B., May, M. J., Pacella, D., Finkenthal, M., Gregory B. C., & Goldstein, W. 
H. 2000, Nuc. Fusion, 40, 847 

French, R. S., Dinerstein, H. L., & Sneden, C. 2000, in "Asymmetrical Planetary Nebulae 
II: From Origins to Microstructures" , ASP Con. 199, eds. J. H. Kastner et al. (San 
Francisco: ASP), 325 

Gallino, R., Arlandini, C, Busso, M., Lugaro, M., Travaglio, C, Straniero, O., Chieffi, A., 
& Limongi, M. 1998, ApJ, 497, 388 

Geballe, T. R., Burton, M. G., & Isaacman, R. 1991, MNRAS, 253, 75 

Goriely, S., & Mowlavi, N. 2000, A&A, 362, 599 

Gould, R. J. 1978, ApJ, 219, 250 

Hanuschik, R. W. 2003, A&A, 407, 1157 

Henry, R. B. C, Kwitter, K. B., & Balick, B. 2004, AJ, 127, 2284 
Herwig, F. 2005, ARA&A, 43, 435 

Hora, J. L., Latter, W. B., & Deutsch, L. K. 1999, ApJS, 124, 195 
Hyung, S., & AUer, L. H. 1997a, MNRAS, 292, 71 
Hyung, S., & AUer, L. H. 1997b, ApJ, 491, 242 



-32- 

Hyung, S., AUer, L. H., & Feibelman, W. A. 1994a, MNRAS, 269, 975 

Hyung, S., AUer, L. H., & Feibelman, W. A. 1994b, PASP, 106, 745 

Hyung, S., AUer, L. H., & Feibelman, W. A. 1997, ApJS, 108, 503 

Hyung, S., AUer, L. H., Feibelman, W. A., & Lee, S.-J. 2001, ApJ, 563, 889 

Hyung, S., Keyes, C. D., & AUer, L. H. 1995, MNRAS, 272, 49 

Iben, I., Jr., & Renzini, A. 1983, ARA&A, 21, 271 

Kallman, T., & Bautista, M. 2001, ApJS, 133, 221 

Kappeler, F., Beer, H., & Wisshak, K. 1989, Rep. Prog. Phys., 52, 945 

Keyes, C. D., AUer, L. H., & Feibelman, W. A. 1990, PASP, 102, 59 

Khromov, G. S. 1989, Space Sci. Rev., 51, 339 

Kingdon, J. B., & Ferland, G. J. 1996, ApJS, 106, 205 

Kramers, H. A. 1923, PhU. Mag., 46, 836 

Kudritzki, R.-R, & Puis, J. 2000, ARA&A, 38, 613 

Lester, D. F., HiU, G. J., Doppman, G., & Froning, C. S. 2000, PASP, 112, 384 
Liu, Y., Liu, X.-W., Luo, S.-G., & Barlow, M. J. 2004a, MNRAS, 353, 1231 
Liu, Y., Liu, X.-W., Barlow, M. J., & Luo, S.-G. 2004b, MNRAS, 353, 1251 
Loch, S. D., et al. 2002, Phys. Rev. A, 66, 052708 

Lugaro, M., Davis, A. M., Gallino, R., Pellin, M. J., Straniero, O., & Kappeler, F. 2003, 
ApJ, 593, 486 

Lumsden, S. L., Puxley, P. J., & Hoare, M. G. 2001, MNRAS, 328, 419 
MazziteUi, G., & Mattioh, M. 2002, Atom. Data Nuc. Data Tab., 82, 313 
Mendez, R. H., Kudritzki, R. P., & Herrero, A. 1992, A&A, 260, 329 

Osterbrock, D. E. 1989, Astrophysics of Gaseous Nebulae and Active Galactic Nuclei (Mill 
Valley: University Science Books) 



-33- 



Osterbrock, D. E., & Ferland, G. J. 2006, Astrophysics of Gaseous Nebulae and Active 
Galactic Nuclei, 2nd Ed. (Sausalito: University Science Books) 

Osterbrock, D. E., Fulbright, J. P., Martel, A. R., Keane, M. J., Trager, S. C., & Basri, G. 
1996, PASP, 108, 277 

Peimbert, M., Peimbert, A., Ruiz, M. T., & Esteban, C. 2004, ApJS, 150, 431 

Pequignot, D., & Baluteau, J.-P. 1988, A&A, 206, 298 

Pequignot, D., & Baluteau, J.-P. 1994, A&A, 283, 593 (PB94) 

Pottasch, S. R., Bernard-Salas, J., Beintema, D. A., & Feibelman, W. A. 2004, A&A, 423, 
593 

Pottasch, S. R., & Surendiranath, R. 2005, A&A, 432, 139 
Ranch, T. 2003, A&A, 403, 709 

Richter, M., Ulm, G., Gerth, C., Tiedtke, K., Feldhaus, J., Sorokin, A. A., Shmaenok, L. 
A., & Bobashev, S. V. 2003, in "X-Ray and Inner Shell Processes", AIP Conf. Proc. 
652, eds. A. Bianconi et al. (New York: AIP), 165 

Rudy, R. J., Lynch, D. K., Mazuk, S., Puetter, R. C., & Dearborn, D. S. P. 2001, AJ, 121, 
362 

Savin, D. W., et al. 1999, ApJS, 123, 687 
Savin, D. W., et al. 2006, ApJ, 642, 1275 
Schoning, T. 1997, A&AS, 122, 277 
Schoning, T., & Butler, K. 1998, A&AS, 128, 581 
Seaton, M. J. 1959, MNRAS, 119, 81 

Sharpee, B., WiUiams, R., Baldwin, J. A., & van Hoof, P. A. M. 2003, ApJS, 149, 157 

Simmerer, J., Sneden, C., Cowan, J. J., Collier, J., Woolf, V. M., & Lawler, J. E. 2004, ApJ, 
617, 1091 

Smits, D. P. 1991, MNRAS, 251, 316 

Sterling, N. C. 2006a, BAAS, 207, #93.03 



-34- 



Sterling, N. C. 2006b, Ph. D. Thesis, University of Texas at Austin 
Sterhng, N. C, & Dinerstein, H. L. 2003a, RMxA&A Ser. de Conf., 18, 133 
Sterhng, N. C, & Dinerstein, H. L. 2003b, BAAS, 203, #11.05 

Sterhng, N. C, & Dinerstein, H. L. 2004, in "Asymmetrical Planetary Nebulae III: Winds, 
Structure and the Thunderbird" , ASP Conf. Ser. 313, eds. M. Meixner et al. 
(San Francisco: ASP), 410 

Sterling, N. C, & Dinerstein, H. L. 2005a, in "Cosmic Abundances as Records of Stellar 
Evolution and Nucleosynthesis", ASP Conf. Ser. 336, eds. T. G. Barnes III & F. N. 
Bash (San Francisco: ASP), 367 

Sterling, N. C, & Dinerstein, H. L. 2005b, RMxA&A Ser. de Conf., 23, 1 

Sterling, N. C, & Dinerstein, H. L. 2006, in "Planetary Nebulae in our Galaxy and Beyond", 
lAU Symp. 234, eds. M. J. Barlow & R. H. Mendez (Cambridge: Cambridge Univ. 
Press), in press (astro-ph/0605059) 

Sterling, N. C, Dinerstein, H. L., & Bowers, C. W. 2002, ApJ, 578, L55 

Sterling, N. C, Dinerstein, H. L., Bowers, C. W., & Redfield, S. 2005, ApJ, 625, 368 

Storey, P. J., & Hummer, D. G. 1995, MNRAS, 272, 41 

Straniero, O., GaUino, R., & Cristallo, S. 2006, Nucl. Phys. A, 777, 311 

Tsamis, Y. G., Barlow, M. J., Liu, X.-W., Danziger, I. J., & Storey, P. J. 2003, MNRAS, 
345, 186 

Travagho, C, GaUino, R., Arnone, E., Cowan, J., Jordan, F., & Sneden, C. 2004, ApJ, 601, 
864 

Treffers, R. R., Fink, U., Larson, H. P., & Gautier, T. N., III. 1976, ApJ, 209, 793 

van Hoof, P. A. M., Weingartner, J. C, Martin, P. G., Volk, K., & Ferland, G. J. 2004, 
MNRAS, 350, 1330 

Voronov, G. S. 1997, Atom. Data Nuc. Data Tab., 65, 1 

Wiese, W. L., Fuhr, J. R., & Deters, T. M. 1996, Atomic Transition Probabilities of Carbon, 
Nitrogen, and Oxygen : A Critical Data Compilation (J. Phys. Chem. Ref. Data, 
Monogr. 7) (Washington, DC: ACS/NIST) 



- 35 - 



Zhang, Y., Liu, X.-W., Luo, S.-G., Pequignot, D., & Barlow, M. J. 2005, A&A, 442, 249 

Zhang, Y., Wilhams, R., Pellegrini, E., Cavagnolo, K., Baldwin, J. A., Sharpee, B., Phillips, 
M., & Liu, X.-W. 2006, in lAU Symp. 234, "Planetary Nebulae in our Galaxy and 
Beyond", eds. M. J. Barlow & R. H. Mendez, in press ( astro-ph/0605181} 



This preprint was prepared with the AAS I^TJrjX macros v5.2. 



-36- 



Table 1. Optical and Near-Infrared Data References 



Object 


Optical 


Spectral 




NIR Data 


Other 


Name 


Data Ref. 


Res. (A)*^ 


c(H/3) 


Ref. 


Ref.'' 


IC418 


Sharpee et al. 2003 


0.1-0.3 


0.34 


Paper II 


Pottasch et al. 2004 


IC 5117 


Hyung ct al. 2001 


0.05 0.15 


1.40 


Rudy ct al. 2001 




NGC6572 


Liu et al. 2004a,b 


1.5-3.5 


0.48 


Paper II 


Hyung et al. 1994a^ 


NGC6741 


Liu et al. 2004a,b 


1.5-3.5 


1.15 


Paper II 




NGC6790 


Liu et al. 2004a,b 


1.5-3.5 


1.10 


Paper II 


AUer et al. 1996^ 


NGC6826 


Liu et al. 2004a,b 


1.5-3.5 


0.06 


Paper II 




NGC6884 


Liu et al. 2004a,b 


1.5-3.5 


1.00 


Paper II 




NGC6886 


Hyung et al. 1995 


0.05-0.15 


0.90 


Paper II 


Pottasch & Surendiranath 2005 


NGC 7027 


Zhang et al. 2005 


1.5-4.5 


1.37 


Geballe et al. 1991 




NGC 7662 


Liu et al. 2004a,b 


1.5-3.5 


0.18 


Geballe et al. 1991 





Note. — Optical and NIR data references are given, along with the optical spectral resolutions and extinction 
coefficients c(H/3). 

^The range of FWHM resolution of the optical spectrum. 

''Reference for abundances, UV, and IR data if not from optical data reference. 

'^[Kr IV] A5867.7 is taken from this reference and scaled to [Kr IV] A5346.0 of Liu et al. (2004a,b). AU other 
abundance, optical, and IR data (except the K band) are from Liu et al. 
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Table 2. Dereddened Se and Kr Line Intensities Relative to H/3 (/(H/5) = 10"^) 



PN 


[Se IV] 


[Kr III] 


[Kr III] 


[Kr IV] 


[Kr IV] 


[Kr V] 


Name 


2.287 iim 


A6826.7'^ 


2.199 /xm 


A5346.0 


A5867.7*^ 


A6256.1 


IC 418 




2.90±0.87'= 


6.99±0.54 


0.35±0.11'= 


0.35±0.ir 




IC 5117 


31.5±3.4 




4.23±1.06 


12.9±3.9 


14.5±4.4 




NGC 6572 


8.87±0.51 


1.80±0.54 


2.49±0.61 


4.70±0.94 


6.40±1.92 




NGC 6741 


9.00±2.23 


3.70±0.37 


5.40±1.38 


6.90±1.38 


5.50±1.10 




NGC 6790 


8.61±0.91 


0.45±0.14 




5.30±1.06 


5.20±1.04 




NGC 6826 


5.41±1.55 


0.93±0.28 




2.30±0.69 


4.50±1.35 




NGC 6884 


20.3±3.0 


0.83±0.25 




9.50±0.95 


11.4±1.1 




NGC 6886 


14.6±1.6 




5.89±1.19 


8.00±2.4^ 


5.40±1.6^ 




NGC 7027 


21.2±1.4 


4.10±1.23'= 


8.74±1.52 


14.8±4.4^ 


22.2±6.7^ 


2.50±0.75^ 


NGC 7662 


23.4±5.2 


0.82±0.25 




15.9±1.6 


16.0±1.6 





Note. — Error bars are 1-a uncertainties, and account only for uncertainties in the 
flux measurements. Line intensities are taken from the references in Table 1. 

^Corrected for He I A6827.9 contamination as described in Appendix A. 3. 

"^Corrected for He II A5869.0 contamination as described in Appendix A. 2. 

*^An uncertainty of 30% in the line intensity is assumed, since errors were not reported 
in the optical data reference. 
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Table 3. Cloudy Input Parameters 



Parameter IC 418^ IC 418^ IC 5117 NGC 6572 NGC 6741 NGC 6790 Solar^ 



Teff (103 K) 


38.8±3.0 


37.3±2.0 


125.25±10.0 


78.71±7.0 


152.49t^°;° 


96.65±10.0 




log g 


3.4 


3.4 


8.0 


5.5 


5.84 


5.5 






4390 


7175 


3300 


3820 


910 


5800 




nn (10^ cm"^) 


13.0 


13.0 


57.4 


17.33 


5.5 


50.75 




Rin (log pc) 


-3.0 


-3.0 


-2.1 


-2.0 


-2.24 


-2.3 




-Rout (log pc) 


-1.34 


-1.45 


-1.61 


-1.44 


-1.26 


-1.74 




D (kpc)'* 


1.6 


2.1 


14.0 


2.6 


6.5 


11.2 




He/H 


0.106 


0.141 


0.100 


0.104 


0.105 


0.094 


0.085 


C/H 


8.78 


8.68 


8.45 


8.75 


8.55 


8.37 


8.39 


N/H 


8.03 


8.03 


7.93 


8.09 


8.32 


7.60 


7.78 


0/H 


8.71 


8.57 


8.40 


8.50 


8.64 


8.38 


8.66 


Ne/H 


8.01 


8.25 


7.70 


7.75 


7.90 


7.61 


7.84 


S/H 


6.60 


6.65 


6.57 


6.36 


6.72 


6.27 


7.14 


Cl/H 


4.93 


4.89 


5.02 


4.89 


4.96 


4.73 


5.50 


Ar/H 


6.21 


6.19 


6.07 


6.09 


6.34 


5.75 


6.18 


Se/H 






3.67±0.05 


3.11±0.08 


3.39±0.12 


3.10±0.05 


3.33 


Kr/H 


3.90±0.21 


3.77±0.23 


3.76±0.11 


3.54±0.18 


3.75±0.13 


3.24±0.14 


3.28 



Note. — Input parameters for final Cloudy models. The stellar effective temperature;, gravity, and 
luminosity (relative to solar), and nebular hydrogen density, inner and outer radii, and implied distance are 
given in the first seven rows. Abundances follow, with He/H on a linear scale, and the other elements given 
as 12+ log(X/H). 

^IC 418 model with Rauch NLTE model atmosphere and no mass loss. 

''IC 418 model with stellar atmosphere including mass-loss (see text, §3.1.1). 

^From Asplund et al. (2005) 

^Distance implied from the model, calculated by comparing the predicted H/3 luminosity to the observed 
H/3 flux from Cahn et al. (1992). Sincx; we did not attempt to optimize L(H/3) in our models, these values 
should not be viewed as accurate estimates of the distances to the modeled PNe. 
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Table 4. Cloudy Input Parameters 



Parameter 


NGC 6826 


NGC 6884 


NGC 6886 


NGC 7027 


NGC 7662 


Solar 


ieflf (^iU J\j 


OO.DitO.U 


11/1 Q94-1 n n 


1 «/i /I -i-on n 

iD4.4ztZU.U 


1 7^? Q+15.0 
1 / 0.0_3o,o 


lot: «4_i n n 

iZO.DltiU.U 




log g 


o.U 


0. / 


7 n 


D.U 


7 

0. / 




T It 


iUDUU 


iozu 


oyu 


4yuu 


0/1 7n 

Z4 / U 




77,H V-'-'-' cm j 


i.D 




S 1 
o.i 


00. u 


O.Zi 




Rin (log Pc) 


— i.D 


O 7 


O 70 


— Z.UD 


— i.OU 




-Rout (log pc) 


— U.oD 


— 1.4D 


— i.Do 


1 R7 
— i.O ( 


— 1.4o 






4 1 


6 1 


5 4 


3 7 


1 7 




Hc/H 


0.099 


0.099 


0.107 


0.100 


0.103 


0.085 


C/H 


8.33 


8.74 


8.89 


9.05 


8.65 


8.39 


N/H 


7.95 


8.33 


8.47 


8.37 


7.97 


7.78 


0/H 


8.71 


8.63 


8.79 


8.82 


8.59 


8.66 


Ne/H 


7.86 


7.84 


8.17 


8.07 


7.81 


7.84 


S/H 


6.33 


6.68 


6.84 


7.00 


6.65 


7.14 


Cl/H 


4.99 


5.02 


5.16 


5.32 


5.14 


5.50 


Ar/H 


6.13 


6.23 


6.43 


6.44 


6.15 


6.18 


Se/H 


2.89±0.13 


3.45±0.10 


3.59±0.11 


3.78±0.12 


3.50±0.10 


3.33 


Kr/H 


3.57±0.20 


3.65±0.16 


3.81±0.18 


4.26±0.12 


3.84±0.14 


3.28 



Note. — Input parameters for final Cloudy models. The stellar effective temperature, 
gravity, and luminosity (relative to solar), and nebular hydrogen density, inner and outer 
radii, and implied distance are given in the first seven rows. Abundances follow, with He/H 
on a linear scale, and the other elements given as 12+ log(X/H). 

^Prom Asplund et al. (2005) 

^Distance implied from the model, calculated by comparing the predicted H/9 luminosity 
to the observed H/3 flux from Cahn et al. (1992). Since we did not attempt to optimize iv(H/3) 
in our models, these values should not be viewed as accurate estimates of the distances to 
the modeled PNe. 
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Table 5. Cloudy Predicted vs. Observed Diagnostics and Line Intensities 



Ion 


X 




IC 418 


a 


IC 418^' 


IC 5117 


NGC 6572 




ir reu. 




Obs. 








WDS. 




ir reu. 


yjos. 








WDS. 


[Ar III]na 




2 


23E+02 


2 


71E+02 


2 


19E+02 


2 


71E+02 


7 


66E+01 


1.35E+02 


1 


03E+02 


1 


69E+02 


[Ar IV] nn 




2 


OOE+00 


1 


20E+00 


2 


03E+00 


1 


20E+00 


4 


lOE+00 


4.18E+00 


2 


23E+00 


1 


98E+00 


[CI III]nn 




5 


07E-01 


5 


lOE-01 


4 


97E-01 


5 


lOE-01 


2 


96E-01 


3.28E-01 


4 


54E-01 


4 


40E-01 


[N II]na 




6 


58E+01 


7 


84E+01 


6 


75E+01 


7 


84E+01 


1 


78E+01 


2.53E+01 


3 


76E+01 


3 


78E+01 


[0 II] na 




4 


99E+00 


6 


98E+00 


4 


98E+00 


6 


98E+00 


8 


71E-01 


6.95E-01 


2 


79E+00 


3 


06E+00 


[O II] im 




3 


26E+00 


2 


37E+00 


3 


28E+00 


2 


37E+00 


3 


64E+00 


3.10E+00 


3 


32E+00 


2 


48E+00 


[O III]na 




3 


12E+02 


3 


08E+02 


2 


98E+02 


3 


08E+02 


6 


83E+01 


8.83E+01 


2 


94E+01': 


4 


30E+01<: 


[S Iljnn 




4 


84E-01 


4 


70E-01 


4 


84E-01 


4 


70E-01 


4 


46E-01 


4.36E-01 


4 


82E-01 


4 


70E-01 


Hi 


4861.3 


1 


OOE+02 


1 


OOE+02 


1 


OOE+02 


1 


OOE+02 


1 


OOE+02 


l.OOE+02 


1 


OOE+02 


1 


OOE+02 


Hi 


2.166/um 


2 


83E+00 


2 


82E+00 


2 


85E+00 


2 


82E+00 


2 


51E+00 


2.64E+00 


2 


64E+00 


2 


71E+00 


He I 


4471.5 


3 


28E+00 


4 


49E+00 


5 


65E+00 


4 


49E+00 


5 


09E+00 


5.30E+00 


5 


52E+00 


5 


83E+00 


He I 


5875.6 


9 


OlE+00 


1 


37E+01 


1 


56E+01 


1 


37E+01 


1 


62E+01 


1.52E+01 


1 


64E+01 


1 


73E+01 


He I 


6678.2 


2 


27E+00 


3 


87E+00 


3 


88E+00 


3 


87E+00 


3 


44E+00 


3.28E+00 


3 


77E+00 


4 


24E+00 


He II 


1640.4 


3 


22E-04 














5 


90E+01 


5.88E+01 


2 


25E+00 


2 


84E+00 


He II 


4685.7 


3 


32E-05 














8 


18E+00 


9.44E+00 


2 


47E-01 


4 


09E-01 


He II 


5411.5 


2 


32E-06 














6 


21E-01 


9.20E-01 


1 


67E-02 


3 


60E-02 


He II 


2.189/im 


2 


57E-07 














5 


97E-02 


7.43E-02 


1 


64E-03 






[C II] 


2325.6 


9 


21E+01 


9 


57E+01 


6 


14E+01 


9 


57E+01 


3 


93E+01 


8.30E+01 


6 


05E+01 


7 


02E+01 


[C III] 


1908.7 


3 


53E+01 


3 


24E+01 


5 


40E+01 


3 


24E+01 


6 


76E+02 


4.18E+02 


5 


60E+02 


3 


lOE+02 


[cIV] 


1549.5 


2 


52E-01 






2 


04E-02 






5 


92E+02 


1.27E+02 


1 


56E+02 






[N I] 


5199.1 


5 


61E-02 


1 


17E-01 


5 


lOE-02 


1 


17E-01 


5 


81E-02 


8.80E-02 


1 


74E-01 


4 


78E-01 


[N II] 


5754.6 


3 


31E+00 


2 


76E+00 


3 


14E+00 


2 


76E+00 


3 


24E+00 


2.46E+00 


2 


65E+00 


2 


62E+00 


[N II] 


6548.0 


5 


51E+01 


5 


36E+01 


5 


36E+01 


5 


36E+01 


1 


46E+01 


1.42E+01 


2 


52E+01 


2 


54E+01 


[N II] 


6583.5 


1 


63E+02 


1 


63E+02 


1 


58E+02 


1 


63E+02 


4 


29E+01 


5.06E+01 


7 


45E+01 


7 


57E+01 


[N III] 


1749.3 


8 


57E-01 


6 


86E-01 


1 


82E+00 


6 


86E-01 


3 


33E+01 


1.24E+01 


2 


41E+01 


1 


54E+01 


[N III] 


57.32/im 


1 


80E+00 


2 


98E+00 


3 


37E+00 


2 


98E+00 


7 


73E-01 




4 


23E+00 


6 


33E+00 


[N IV] 


1484.9 


1 


32E-04 






1 


06E-07 






3 


53E+01 




6 


69E+00 






[nV] 


1240.8 


















1 


62E+00 




6 


18E-03 






[Ol] 


6300.3 


1 


85E+00 


2 


18E+00 


1 


15E+00 


2 


18E+00 


4 


26E+00 


6.94E+00 


4 


53E+00 


7 


09E+00 


[O I] 


6363.8 


5 


91E-01 


7 


59E-01 


3 


66E-01 


7 


59E-01 


1 


36E+00 


2.29E+00 


1 


44E+00 


2 


33E+00 


[Oil] 


3726.0 


1 


84E+02 


1 


24E+02 


1 


18E+02 


1 


24E+02 


1 


31E+01 


8.67E+00 


3 


56E+01 


2 


89E+01 


[Oil] 


3728.8 


5 


66E+01 


5 


23E+01 


3 


60E+01 


5 


23E+01 


3 


59E+00 


2.80E+00 


1 


07E+01 


1 


16E+01 


[Oil] 


7319.5 


2 


75E+01 


1 


39E+01 


1 


76E+01 


1 


39E+01 


1 


09E+01 


8.80E+00 


9 


47E+00 


8 


30E+00 


[Oil] 


7330.2 


2 


08E+01 


1 


15E+01 


1 


33E+01 


1 


15E+01 


8 


21E+00 


7.69E+00 


7 


16E+00 


6 


68E+00 


[oIII] 


1661.6 


7 


73E-01 






8 


37E-01 






5 


36E+01 


3.80E+01 


2 


43E+01 


1 


19E+01 


[oIII] 


4363.2 


9 


37E-01 


9 


35E-01 


9 


81E-01 


9 


35E-01 


2 


83E+01 


1.97E+01 


1 


52E+01 


9 


14E+00 


[oIII] 


4958.9 


7 


28E+01 


7 


27E+01 


7 


28E+01 


7 


27E+01 


4 


82E+02 


3.95E+02 


4 


46E+02 


4 


06E+02 


[oIII] 


5006.8 


2 


19E+02 


2 


15E+02 


2 


19E+02 


2 


15E+02 


1 


45E+03 


1.35E+03 


1 


34E+03 






[oIII] 


51. 80^111 


1 


05E+01 


1 


37E+01 


9 


81E+00 


1 


37E+01 


6 


93E+00 




2 


45E+01 


4 


06E+01 


[oIV] 


1400.7 


3 


22E-07 














5 


02E+00 




9 


21E-03 






[oIV] 


25. 88^111 


2 


27E-05 














5 


07E+00 




2 


16E-01 






[Ne II] 


12.81^tm 


4 


70E+01 


4 


77E+01 


1 


lOE+02 


4 


77E+01 


5 


16E-01 




1 


05E+00 






[Ne III] 


3869.1 


3 


49E+00 


3 


09E+00 


1 


67E+00 


3 


09E+00 


1 


54E+02 


1.04E+02 


1 


14E+02 


9 


99E+01 


[Ne III] 


3967.8 


1 


05E+00 


9 


75E-01 


5 


02E-01 


9 


75E-01 


4 


65E+01 


4.26E+01 


3 


44E+01 


3 


05E+01 
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Table 5 — Continued 



Ion 


A 


IC 418'=- 


IC 418*= 


IC 5117 


NGC 6572 




Pred. 




Obs. 




Pred. 


Obs. 




Pred. 




Obs. 




Pred. 




Obs. 


[Ne III] 


15.55/im 


6 


83E+00 


8 


54E+00 


3 


03E+00 


8.54E+00 


5 


56E+01 






7 


88E+01 


4 


82E+00 


[Ne IV] 


2423.8 


7 


78E-09 












8 


96E+00 


1 


llE+01 


1 


13E-02 






[Nc IV] 


4720.0 


2 


49E-11 












2 


31E-01 


1 


lOE-01 


9 


25E-05 






[Nc V] 


3425.5 
















4 


09E-01 






5 


llE-07 






[Nc V] 


14.32^111 
















2 


95E-01 






1 


81E-06 






[Ne V] 


24.31/im 
















5 


17E-02 






5 


06E-07 






[S II] 


4068.6 


2 


51E+00 


1 


79E+00 


2 


16E+00 


1.79E+00 


4 


79E+00 


3 


03E+00 


2 


20E+00 


1 


97E+00 


[S II] 


4076.3 


8 


15E-01 


7 


65E^01 


6 


99E-01 


7.65E-01 


1 


53E+00 


1 


13E+00 


7 


13E-01 


6 


31E-01 


[S II] 


6716.4 


1 


86E+00 


2 


08E+00 


1 


60E+00 


2.08E+00 


8 


06E-01 


9 


OOE-01 


1 


16E+00 


1 


OlE+00 


[S II] 


6730.8 


3 


85E+00 


4 


42E+00 


3 


31E+00 


4.42E+00 


1 


81E+00 


2 


06E+00 


2 


41E+00 


2 


13E+00 


[S III] 


6312.1 


1 


17E+00 


8 


57E-01 


1 


25E+00 


8.57E-01 


2 


07E+00 


2 


67E+00 


1 


OlE+00 


1 


29E+00 


[S III] 


9068.6 


2 


04E+01 


1 


78E+01 


2 


26E+01 


1.78E+01 


1 


43E+01 


2 


23E+01 


9 


91E+00 






[S III] 


9530.6 


5 


05E+01 


4 


23E+01 


5 


60E+01 


4.23E+01 


3 


55E+01 


6 


57E+01 


2 


46E+01 






[S III] 


18.71/im 


1 


56E+01 


1 


36E+01 


1 


70E+01 


1.36E+01 


2 


42E+00 






4 


61E+00 


1 


43E+00 


[s IV] 


10.5lAtm 


1 


14E+00 


1 


14E+00 


1 


70E+00 


1.14E+00 


1 


95E+01 


3 


36E+01 


2 


17E+01 






[CI II] 


8578.7 


3 


35E-01 


2 


84E-01 


2 


65E-01 


2.84E-01 


2 


63E-01 






2 


12E-01 






[CI III] 


5517.7 


1 


88E-01 


1 


82E-01 


1 


84E-01 


1.82E-01 


1 


61E-01 


1 


40E-01 


2 


36E-01 


2 


61E-01 


[CI III] 


5537.9 


3 


70E-01 


3 


56E-01 


3 


70E-01 


3.56E-01 


5 


42E-01 


4 


30E-01 


5 


20E-01 


5 


93E-01 


[CI IV] 


7530.5 


2 


18E-01 






2 


20E-01 




3 


89E-01 


4 


20E-01 


2 


89E-01 


1 


17E-01 


[CI IV] 


8045.6 


4 


38E-01 






4 


41E-01 




7 


80E-01 


9 


40E-01 


5 


78E-01 






[Ar II] 


6. 99/1111 


7 


61E+00 


5 


39E+00 


4 


71E+00 


5.39E+00 


2 


72E-01 






2 


60E-01 






[Ar III] 


5191.8 


4 


62E-02 


3 


86E-02 


5 


75E-02 


3.86E-02 


2 


03E-01 


1 


60E-01 


1 


55E-01 


1 


21E-01 


[Arlll] 


7135.8 


8 


32E+00 


8 


26E+00 


1 


02E+01 


8.26E+00 


1 


25E+01 


1 


68E+01 


1 


28E+01 


1 


78E+01 


[Arlll] 


7751.1 


2 


OlE+00 


2 


20E+00 


2 


45E+00 


2.20E+00 


3 


03E+00 


4 


28E+00 


3 


lOE+00 


4 


18E+00 


[Ar III] 


8.99^111 


7 


85E+00 


6 


59E+00 


9 


42E+00 


6.59E+00 


4 


92E+00 






7 


25E+00 






[Ar i V ] 


4711.4 


4 


21E-03 


3 


OOE-03 


9 


56E-04 


3.00E-03 


1 


58E+00 


1 


13E+00 


1 


36E+00 


1 


lOE+00 


[Ar IV] 


4740.2 


8 


43E-03 


3 


60E-03 


1 


94E-03 


3.60E-03 


6 


49E+00 


4 


72E+00 


3 


02E+00 


2 


18E+00 


[Ar IV] 


7170.6 


1 


17E-04 






2 


81E-05 




2 


59E-01 


2 


90E-01 


7 


19E-02 


5 


lOE-02 


[Ar V] 


6435.1 


9 


88E-10 












1 


18E-01 


3 


OOE-02 


2 


16E-04 






[Ar V] 


7005.7 


2 


12E-09 












2 


54E-01 


3 


OOE-01 


4 


62E-04 






[Ar V] 


7.90Mm 


1 


08E-08 












3 


49E-01 






1 


26E-03 






[Ar VI] 


4.54/im 
















1 


78E-01 






4 


14E-06 






[Se IV] 


2.287/im 


2 


67E-02 






2 


23E-01 




3 


15E-01 


3 


15E-01 


8 


87E-02 


8 


87E-02 


[Kr III] 


6826.7 


3 


30E-02 


2 


90E-02 


3 


26E-02 


2.90E-02 


2 


14E-02 






1 


29E-02 


1 


80E-02 


[Kr III] 


2.199/im 


1 


07E-01 


6 


99E-02 


9 


50E-02 


6.99E-02 


3 


54E-02 


4 


23E-02 


2 


48E-02 


2 


49E-02 


[Kr IV] 


5346.0 


2 


OlE-03 


3 


50E-03 


1 


77E-03 


3.50E-03 


1 


26E-01 


1 


29E-01 


4 


94E-02 


4 


70E-02 


[Kr IV] 


5867.7 


3 


02E-03 


3 


50E-03 


2 


62E-03 


3.50E-03 


1 


82E-01 


1 


45E-01 


7 


04E-02 


6 


40E-02 


[Kr V] 


62-56.1 


3 


41E-09 












2 


35E-03 






1 


12E-04 






[Kr V] 


821:11 




12E-()!) 












2 








1 


()2E-()i 







Note. — Comparison of predicted and observed line intensities relative to H,(3 (on the scale /(Hfl) = 100) and diagnostic 
ratios used in the spectral fits. Wavelengths are in air for A > 3000 A and vacuum otherwise. Listed wavelengths are in A 
except where noted. The diagnostic ratios are defined as follows: [Ar III]na = [Ar III] (7135.8 + 7751.1)/5191.8 ; [Ar IV]nn 
= [Ar IV] 4740.2/4711.4; [CI III]nn = [CI III] 5517.7/5537.9; [N II]iia = [N II] (6548.0 + 6583.5)/5754.6; [O II]na = [O II] 
(3726.0 + 3728.8)/(7319.5 + 7330.2); [O II]im = [O II] 3726.0/3728.8; [O III]na = [O III] (4958.9 + 5006.8)/4363.2; and 
[S Iljnn = [S II] 6716.4/6730.8. 
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*IC 418 model with Rauch NLTE model atmosphere and no mass loss. 

''IC 418 model with stellar atmosphere including mass-loss (see text, §3.1.1). 

■=[0 III]na = [O III] 4958.9/4363.2. 
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Table 6. Cloudy Predicted vs. Observed Diagnostics and Line Intensities 



Ion 



NGC 6741 



Pred. 



Obs. 



NGC 6790 



Pred. 



Obs. 



NGC 6826 



Prod. 



Obs. 



NGC 6884 



Pred. 



Obs. 



[Ar III]na 
[Ar IV] nn 
[CI III]nn 
[N II] na 
[O II]iia 
[O Ilinn 
[O Illjna 
[S Iljnn 
Hi 
Hi 
He I 
He I 
He I 
He II 
He II 
He II 
He II 
[ClI] 
[ClII] 

[CIV] 

[N I] 
[N II] 
[N II] 
[N II] 
[N III] 
[N III] 
[N IV] 

[nV] 

[O I] 

[O I] 

[Oil] 

[Oil] 

[Oil] 

[Oil] 

[OIII] 

[oIII] 

[oIII] 
[oIII] 
[oIII] 
[oIV] 
[OIV] 

[Nc II] 
[No III] 
[Ne III] 
[Ne III] 



4861.3 

2.166^tm 
4471.5 
5875.6 
6678.2 
1640.4 
4685.7 
5411.5 

2.189/^m 
2325.6 
1908.7 
1549.5 
5199.1 
5754.6 
6548.0 
6583.5 
1749.3 

57.32^m 
1484.9 
1240.8 
6300.3 
6363.8 
3726.0 
3728.8 
7319.5 
7330.2 
1661.6 
4363.2 
4958.9 
5006.8 

51. 80^111 
1400.7 

25.88/um 

12.81;um 
3869.1 
3967.8 

15.55/im 



9.64E+01 

1.28E+00 

7.47E-01 

5.66E+01 

7.58E+00 

2.63E+00 

2.39E+01'' 

5.81E-01 

l.OOE+02 

2.56E+00 

3.62E+00 

1.05E+01 

2.49E+00 

2.88E+02 

4.06E+01 

3.09E+00 

2.99E-01 

6.87E+01 

4.72E+02 

5.31E+02 

4.23E+00 

7.08E+00 

l.OlE+02 

2.99E+02 

5.73E+01 

1.78E+01 

7.41E+01 

5.98E+01 

1.59E+01 

5.07E+00 

1.58E+02 

5.99E+01 

1.63E+01 

1.24E+01 

4.09E+01 

2.08E+01 

4.97E+02 

1.50E+03 

6.03E+01 

2.54E+01 

2.05E+02 

2.44E+00 

1.50E+02 

4.51E+01 

9.76E+01 



1.32E+02 

1.26E+00 

7.30E-01 

6.48E+01 

9.47E+00 

2.04E+00 

2.67E+01'' 

5.70E-01 

l.OOE+02 

2.61E+00 

4.58E+00 

1.30E+01 

3.18E+00 

2.41E+02 

3.62E+01 

2.79E+00 

2.75E-01 

1.24E+02 

5.10E+02 

3.25E+02 

2.71E+00 

6.39E+00 

1.02E+02 

3.12E+02 

5.84E+01 

1.08E+01 

2.06E+01 

7.04E+01 

1.66E+01 

5.41E+00 

1.51E+02 

7..39E+01 

1.31E+01 

1.06E+01 

3.36E+01 

1.89E+01 

5.05E+02 

4.03E+01 
5.50E+01 
1.79E+02 
8.29E+00 
1.64E+02 
3.82E+01 
1.25E+02 



7.52E+01 
3.85E+00 
3.07E-01 
1.90E+01 
9.64E-01 
3.62E+00 
1.77E+01 
4.48E-01 
l.OOE+02 
2.54E+00 
5.08E+00 
1.61E+01 
3.40E+00 
1.99E+01 
2.64E+00 
1.99E-01 
1.90E-02 
2.72E+01 
4.04E+02 
4.70E+02 
1.51E-02 
1.14E+00 
5.47E+00 
1.61E+01 
1.52E+01 
3.90E-01 
1.52E+01 
3.42E-01 
3.14E+00 
l.OOE+00 
1.13E+01 
3.12E+00 
8.53E+00 
6.42E+00 
5.33E+01 
2.78E+01 
4.91E+02 
1.48E+03 
7.89E+00 
2.27E+00 
2.28E+00 
3.71E-01 
1.30E+02 
3.90E+01 
4.85E+01 



9.16E+01 
4.06E+00 

4.07E-01 
2.27E+01 
1.51E+00 
2.45E+00 
2.24E+01 
4.57E-01 
l.OOE+02 
2.58E+00 
6.42E+00 
1.71E+01 
4.07E+00 
2.30E+01 
3.83E+00 
2.80E-01 
3.43E-02 
5.28E+01 
3.36E+02 
3.26E+00 
1.66E-01 
1.03E+00 
5.90E+00 
1.76E+01 
1.41E+01 



2.92E+01 
4.52E+00 
1.56E+00 
1.23E+01 
5.01E+00 
6.27E+00 
5.25E+00 
3.44E+01 
2.15E+01 
4.81E+02 

7.57E+00 

1.49E+00 

1.27E+02 
3.70E+01 
4.31E+01 



2.11E+02 

9.04E-01 

1.08E+00 

9.01E+01 

2.32E+01 

1.85E+00 

7.14E+01=' 

7.21E-01 

l.OOE+02 

2.81E+00 

4.83E+00 

1.27E+01 

3.35E+00 

1.66E-01 

1.17E-02 

6.52E-04 

6.24E-05 

2.55E+00 

3.85E+01 

2.60E+00 

4.42E-03 

1.41E-01 

3.21E+00 

9.47E+00 

2.71E+00 

2.87E+01 

1.72E-01 

5.62E-05 

1.77E-02 

5.66E-03 

1.92E+01 

1.04E+01 

7.17E-01 

5.58E-01 

4.70E+00 

4.92E+00 

3.51E+02 

1.06E+03 

2.02E+02 

5.03E-05 

3.39E-02 

2.30E+00 

5.39E+01 

1.63E+01 

l.OlE+02 



2.56E+02 

9.32E-01 

l.lOE+00 

7.91E+01 

1.48E+01 

1.48E+00 

6.71E+01'' 

7.30E-01 

l.OOE+02 

2.82E+00 

5.19E+00 

1.43E+01 

4.07E+00 

2.70E-02 



5.65E+00 
2.44E+01 



1.50E-01 
3.04E+00 
8.79E+00 
3.17E+00 
3.13E+01 



6.70E-02 
2.40E-02 
1.49E+01 

l.OlE+01 
9.21E-01 
7.65E-01 
3.33E+00 
3.74E+00 
2.51E+02 

1.68E+02 



4.92E+01 
1.43E+01 
8.41E+01 



l.llE+02 

1.42E+00 

6.70E-01 

4.75E+01 

5.83E+00 

2.89E+00 

3.20E+01'' 

5.21E-01 

l.OOE+02 

2.66E+00 

4.33E+00 

1.23E+01 

3.01E+00 

9.99E+01 

1.42E+01 

1.08E+00 

1.06E-01 

2.24E+01 

4.12E+02 

3.62E+02 

1.37E-02 

8.68E-01 

1.04E+01 

3.08E+01 

3.23E+01 

1.54E+01 

3.64E+01 

1.16E+01 

l.OlE-01 

3.23E-02 

2.10E+01 

7.28E+00 

2.75E+00 

2.10E+00 

2.88E+01 

1.80E+01 

5.76E+02 

1.73E+03 

7.10E+01 

7.94E+00 

7.83E+01 

5.36E-01 

1.25E+02 

3.77E+01 

9.94E+01 



1.51E+02 

1.56E+00 

7.04E-01 

5.25E+01 

7.02E+00 

2.17E+00 

3.78E+01'' 

5.32E-01 

l.OOE+02 

2.70E+00 

4.99E+00 

1.36E+01 

3.08E+00 

l.llE+02 

1.83E+01 

1.33E+00 

1.80E-01 

2.67E+01 

2.43E+02 

4.29E+01 

2.75E-01 

9.14E-01 

1.19E+01 

3.60E+01 

1.76E+01 

1.16E+01 



2.31E+00 
7.72E-01 
2.04E+01 

9.43E+00 
2.27E+00 
2.00E+00 
2.12E+01 
1.25E+01 
4.73E+02 

7.34E+01 

7.70E+01 

1.27E+02 
3.93E+01 
8.28E+01 
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Table 6 — Continued 



Ion 


\ 


NGC 6741 


NGC 6790 


NGC 6826 


NGC 6884 




Pred. 


Obs. 




Pred. 




Obs. 




Pred. 




Obs. 




Pred. 




Obs. 


[Nc IV] 


2423.8 


5 


62E+01 


1.23E+02 


2 


80E+00 






4 


56E-05 








2 


02E+01 


2 


54E4-01 


[Nc IV] 


4720.0 


6 


22E-01 




7 


lOE-02 






1 


42E-07 








2 


24E-01 


1 


89E-01 


[Nc V] 


3425.5 


4 


63E+01 




1 


72E-02 
















1 


26E+00 






[Ne V] 


14.32/im 


9 


06E+01 


3.37E+01 


1 


24E— 02 
















2 


30E+00 






[NeV] 


24.31/nm 


4 


98E+01 


3.82E+01 


2 


32E-03 
















1 


15E4-00 






[S II] 


4068.6 


8 


07E+00 


8.80E+00 


1 


89E+00 


1 


66E+00 


1 


33E— 01 


1 


41E- 


01 


1 


32E+00 


1 


48E+00 


[S II] 


4076.3 


2 


61E+00 


2.89E+00 


6 


07E-01 


5 


06E-01 


4 


29E-02 


1 


75E- 


01 


4 


29E-01 


4 


48E-01 


[S II] 


6716.4 


1 


23E+01 


1.40E+01 


3 


5gE— 01 


3 


97E-01 


5 


19E-01 


4 


26E- 


01 


1 


33E4-00 


1 


31E4-00 


[S II] 


6730.8 


2 


llE+01 


2.43E+01 


7 


95E-01 


8 


71E-01 


7 


20E-01 


5 


82E- 


01 


2 


55E+00 


2 


45E+00 


[S III] 


6312.1 


2 


73E+00 


3.13E+00 


9 


55E— 01 


1 


lOE+00 


4 


41E— 01 


6 


04E- 


01 


1 


72E+00 


1 


36E+00 


[S III] 


9068.6 


2 


53E+01 


2.55E+01 


6 


61E+00 


8 


24E+00 


8 


27E+00 








1 


85E+01 






[S III] 


9530.6 


6 


27E+01 


7.33E+01 


1 


64E+01 


2 


13E+01 


2 


05E+01 








4 


59E+01 






[S III] 


18.7lMm 


2 


15E+01 


2.22E+01 


1 


23E+00 






1 


46E+01 


1 


25E4 


01 


1 


52E4-01 


1 


43E4-01 


[S IV] 


10.51/im 


5 


58E+01 


2.38E+01 


1 


15E+01 


1 


02E+01 


2 


79E+01 


2 


62E4 


01 


7 


12E4-01 


9 


58E4-01 


[CI II] 


8578.7 


4 


84E-01 


5.44E-01 


1 


08E-01 






6 


64E-02 








1 


04E-01 






[CI III] 


5517.7 


4 


73E-01 


6.92E-01 


8 


36E-02 


9 


50E-02 


4 


22E-01 


4 


29E- 


01 


3 


90E-01 


4 


45E-01 


[CI III] 


5537.9 


6 


33E-01 


9.49E-01 


2 


72E-01 


2 


35E-01 


3 


92E-01 


3 


91E- 


01 


5 


82E-01 


6 


34E-01 


[CI IV] 


7530.5 


3 


41E-01 


1.79E-01 


1 


86E-01 


1 


86E— 01 


2 


56E-01 


4 


lOE- 


02 


3 


17E— 01 


2 


89E— 01 


[ciIV] 


8045.6 


6 


82E-01 


4.86E-01 


3 


73E-01 


4 


OOE-01 


5 


14E-01 








6 


35E-01 






[Ar II] 


6. 99^111 


1 


20E+00 




9 


57E-02 






8 


63E-02 








1 


02E— 01 






[Ar III] 


5191.8 


3 


06E— 01 


2.61E— 01 


9 


02E— 02 


9 


OOE— 02 


5 


73E— 02 


5 


60E- 


02 


1 


57E— 01 


1 


20E— 01 


[Ar III] 


7135.8 


2 


38E+01 


2.81E+01 


5 


46E+00 


6 


60E+00 


9 


74E+00 


1 


15E+01 


1 


40E4-01 


1 


54E+01 


[Ar III] 


7751.1 


5 


74E+00 


6.45E+00 


1 


32E+00 


1 


61E+00 


2 


35E+00 


2 


70E4 


00 


3 


39E4-00 


3 


39E+00 


[Ar III] 


8.99Atin 


1 


40E+01 


9.41E+00 


2 


18E+00 


1 


57E+00 


9 


64E+00 


6 


24E4 


00 


8 


85E4-00 


1 


14E4-01 


[Ar IV] 


4711.4 


4 


51E+00 


4.39E+00 


8 


67E-01 


7 


25E-01 


7 


77E-01 


3 


85E- 


01 


3 


34E4-00 


3 


62E+00 


[Ar IV] 


4740.2 


5 


79E+00 


5.86E+00 


3 


34E+00 


2 


94E+00 


7 


02E-01 


3 


59E- 


01 


4 


74E+00 


5 


67E+00 


[Ar IV] 


7170.6 


1 


86E-01 


3.22E-01 


1 


24E-01 


9 


70E-02 


1 


06E-02 








1 


17E-01 


1 


47E-01 


[Ar V] 


6435.1 


1 


llE+00 


5.85E-01 


2 


34E-02 






2 


52E-06 








3 


22E-01 


1 


07E-01 


[Ar V] 


7005.7 


2 


39E+00 


1.27E+00 


5 


02E-02 






5 


40E-06 








6 


90E-01 


1 


94E-01 


[Ar V] 


7.90/im 


4 


33E+00 




6 


84E-02 






2 


27E-05 








1 


37E4-00 






[Ar VI] 


4. 54^111 


9 


93E+00 


1.85E+00 


2 


74E-02 






2 


24E-12 








1 


95E4-00 






[Sc IV] 


2.287^tm 


9 


OOE-02 


9.00E-02 


8 


61E-02 


8 


61E-02 


5 


41E-02 


5 


41E- 


02 


2 


03E-01 


2 


03E-01 


[Kr III] 


6826.7 


2 


42E-02 


3.70E-02 


5 


58E-03 


4 


50E-03 


9 


89E-03 


9 


30E- 


03 


1 


32E-02 


8 


30E-03 


[Kr III] 


2.199//m 


4 


70E-02 


5.40E-02 


9 


15E-03 






2 


82E-02 








2 


69E-02 






[Kr IV] 


5346.0 


5 


96E-02 


6.90E-02 


4 


02E-02 


5 


30E-02 


2 


61E-02 


2 


30E- 


02 


7 


12E-02 


9 


50E-02 


[Kr IV] 


5867.7 


8 


26E-02 


5.50E-02 


5 


78E-02 


5 


20E-02 


3 


80E-02 


4 


50E- 


02 


1 


OOE-01 


1 


14E-01 


[Kr V] 


6256.1 


9 


83E-03 




6 


17E-04 






3 


20E-06 








2 


17E-03 






[Kr V] 


8243.4 


8 


98E-03 


5.50E-02 


5 


64E-04 






2 


92E-06 








1 


98E-03 







Note. — Comparison of predicted and observed line intensities relative to H/3 (on the scale /(H/3) = 100) and diagnostic 
ratios used in the spectral fits. Wavelengths are in air for A > 3000 A and vacuum otherwise. Listed wavelengths are in A 
except where noted. The diagnostic ratios are defined as follows: [Ar III]na = [Ar III] (7135.8 + 7751.1)/5191.8 ; [Ar IV]nn 
= [Ar IV] 4740.2/4711.4; [CI III]nn = [CI III] 5517.7/5537.9; [N II]na = [N 11] (6548.0 + 6583.5)/5754.6; [O II]na = [O II] 
(3726.0 4- 3728.8)/(7319.5 4- 7330.2); [O Iljnn = [O II] 3726.0/3728.8; [O III]na = [O III] (4958.9 4- 5006.8)/4363.2; and 
[S II]nn = [S II] 6716.4/6730.8. 

^[O III]na = [O III] 4958.9/4363.2. 
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Table 7. Cloudy Predicted vs. Observed Diagnostics and Line Intensities 
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Table 7— Continued 
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70E-01 


8 


20E-01 
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1 
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5 


98E+00 


8 


29E+00 


1 


20E+01 
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5 
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3 
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5 
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Tabic 7- 


-Contimicd 








Ion 


A 


NGC 6886 


NGC 7027 


NGC 7662 


Pred. 


Obs. 


Pred. 


Obs. 


Pred. 


Obs. 


[ArV] 


6435.1 


1.18E+00 


l.llE+00 


2.01E+00 


1.31E+00 


6.99E-01 


4.15E-01 


[Ar V] 


7005.7 


2.53E+00 


1.89E+00 


4.31E+00 


3.09E+00 


1.50E+00 


8.53E-01 


[Ar V] 


7.90/ini 


5.32E+00 




6.58E+00 


4.67E+00 


2.90E+00 


2.11E+00 


[Ar VI] 


4.54/xm 


1.29E+01 


7.95E+00 


2.59E+01 


1.33E+01 


3.12E+00 


1.46E+00 


[Se IV] 


2.287/im 


1.46E-01 


1.46E-01 


2.11E-01 


2.12E-01 


2.34E-01 


2.34E-01 


[Kr III] 


6826.7 


2.43E-02 




4.27E-02 


4.10E-02 


1.03E-02 


8.20E-03 


[Kr III] 


2.199/xm 


5.24E-02 


5.89E-02 


8.60E-02 


8.74E-02 


2.12E-02 




[Kr IV] 


5346.0 


5.68E-02 


8.00E-02 


1.46E-01 


1.48E-01 


1.30E-01 


1.59E-01 


[Kr IV] 


5867.7 


7.99E-02 


5.40E-02 


2.14E-01 


2.22E-01 


1.82E-01 


1.60E-01 


[Kr V] 


6256.1 


9.35E-03 




2.80E-02 


2.50E-02 


l.OlE-02 




[KrV] 


8243.4 


8.54E-03 




2.56E-02 




9.21E-03 





Note. — Comparison of predicted and observed line intensities relative to H/3 (on the scale 
/(H/3) = 100) and diagnostic ratios used in the spectral fits. Wavelengths arc in air for A > 3000 A 
and vacuum otherwise. Listed wavelengths are in A except where noted. The diagnostic ratios 
are defined as follows: [Ar III]na = [Ar III] (7135.8 + 7751.1)/5191.8 ; [Ar IV]nn = [Ar IV] 
4740.2/4711.4; [CI III]nn = [CI III] 5517.7/5537.9; [N II]na = [N II] (6548.0 + 6583.5)/5754.6; 
[O II]na= [Oil] (3726.0 + 3728.8)/(7319.5+ 7330.2); [O II]nn = [O II] 3726.0/3728.8; [O III]na 
= [O III] (4958.9 + 5006.8)/4363.2; and [S II]nn = [S II] 6716.4/6730.8. 

=^[0 III]na = [O III] 4958.9/4363.2. 
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Table 8. Se and Kr Abundance Uncertainties from Various Sources 



PN 12+ Se Abundance Errors 12+ Kr Abundance Errors 



Name 


Log(Se/H) 


/ 






Total Err. 


Log(Kr/H) 


/ 




riH 


Total Err. 


IC 418^ 












3.90 


0.144 


0.17 


0.03 


0.21 


IC 418'' 












3.77 


0.179 


0.14 


0.10 


0.23 


IC 5117 


3.67 


0.044 


0.02 


0.02 


0.05 


3.76 


0.111 


0.01 


0.01 


0.11 


NGC 6572 


3.11 


0.025 


0.01 


0.077 


0.08 


3.54 


0.175 


0.01 


0.06 


0.18 


NGC 6741 


3.39 


0.096 


0.079 


0.01 


0.12 


3.75 


0.129 


0.02 


0.01 


0.13 


NGC 6790 


3.10 


0.044 


0.015 


0.015 


0.05 


3.24 


0.134 


0.01 


0.05 


0.14 


NGC 6826 


2.89 


0.109 


0.06 


0.05 


0.13 


3.57 


0.134 


0.109 


0.123 


0.20 


NGC 6884 


3.45 


0.06 


0.015 


0.086 


0.10 


3.65 


0.122 


0.01 


0.12 


0.16 


NGC 6886 


3.59 


0.045 


0.07 


0.079 


0.11 


3.81 


0.16 


0.088 


0.05 


0.18 


NGC 7027 


3.78 


0.028 


0.119 


0.01 


0.12 


4.26 


0.113 


0.057 


0.02 


0.12 


NGC 7662 


3.50 


0.087 


0.061 


0.01 


0.10 


3.84 


0.134 


0.02 


0.046 


0.14 



Note. — 1-a abundance uncertainties from various sources are given, along with the total error estimate. All 
uncertainties are in dex. 

"^IC 418 model with Rauch NLTE model atmosphere and no mass loss. 

''IC 418 model with stellar atmosphere with mass-loss (see text, §3.1.1). 



Table 9. Line Intensity Errors from Atomic Data Uncertainties 



PN /([SelV]) /([KrIIl]) ^([Kr IV]) /([KrV]) 

Name 2.287 //m AA6826.7, 21990 AA5346.0, 5867.7 A6256.1 



CT Rate Coefficient Varied 

IC 418 0.047 0.031 0.044 ^ 

NGC 6741 0.106 0.125 0.075 0.057 

NGC 6884 0.074 0.152 0.052 0.036 

NGC 7027 0.095 0.147 0.071 0.026 

Low-T DR Rate Coefficient Varied 

IC 418 0.196 0.025 0.204 ^ 

NGC 6741 0.095 0.061 0.070 0.114 

NGC 6884 0.081 0.109 0.053 0.120 

NGC 7027 0.080 0.088 0.066 0.137 

PI Cross-Section Varied 

IC 418 0.215 0.090 0.248 ^ 

NGC 6741 0.126 0.136 0.091 0.120 

NGC 6884 0.132 0.207 0.089 0.217 

NGC 7027 0.122 0.177 0.096 0.156 

RR Rate Coefficient Varied 

IC 418 0.067 0.045 0.120 ^ 

NGC 6741 0.030 0.021 0.040 0.089 

NGC 6884 0.023 0.058 0.018 0.102 

NGC 7027 0.034 0.032 0.050 0.118 

Cumulative Uncertainties from Atomic Data Errors 

IC 418 0.274 0.104 0.308 ^ 

NGC 6741 0.177 0.182 0.134 0.180 

NGC 6884 0.162 0.257 0.112 0.249 

NGC 7027 0.165 0.227 0.136 0.218 



Note. — Predicted line intensity errors (in dex) from atomic data un- 
certainties. See ^j3.1.3. 
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Table 10. Values of Physical Parameters in Grid of Cloudy models 



Teff (103 K) 


50 


55 


60 


65 70 75 


80 


85 


90 


95 


100 


110 


120 


130 140 150 


160 


170 


180 




//ii (10''' ciir''') 


1 


2 


3 


4 5 6 


7 


8 


9 


10 


12.5 


15 


20 


25 30 40 


50 


60 


70 


80 








90 100 










\og{U) 


-3.0 


-2.0 


-1.5 


-1.0 -0.5 0.0 


0.5 


1.0 


2.0 





Note. — Values of the input parameters varied in our grid of 
Cloudy models are hsted. All combinations of the listed param- 
eters were input into our grid, for a total of 3761 models. 
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Table 11. Comparison of ICF- and Model-Derived Kr Abundances 



PN 


n(Kr++)/ 




12+ Log(Kr/H) 


12+ Log(Kr/H) 


Name 


n(H+)xl09 


ICF(Kr)^ 


ICF 


Model 


IC 418 


3.18±0.79 


2.05±1.36 


3.81±0.33 


3.90+0.21 


IC 5117 


1.36±0.37 


4.75±1.96 


3.81±0.24 


3.76+0.11 


NGC 6572 


1.16±0.30 


3.10±2.13 


3.56±0.40 


3.54+0.18 


NGC 6741 


1.80±0.29 


4.74±2.42 


3.93±0.26 


3.75+0.13 


NGC 6790 


0.24±0.13 


5.08±2.45 


3.08±0.40 


3.24+0.14 


NGC 6826 


1.03±0.44 


2.16±1.79 


3.35±0.48 


3.57+0.20 


NGC 6884 


0.62±0.23 


5.35±2.32 


3.52±0.30 


3.65+0.16 


NGC 6886 


1.67±0.37 


5.33±2.51 


3.95±0.25 


3.81+0.18 


NGC 7027 


2.42±0.46 


6.13±2.69 


4.17±0.23 


4.26+0.12 


NGC 7662 


0.40±0.14 


7.98±3.14 


3.50±0.25 


3.84+0.14 



Note. — Comparison of Kr abundances derived from ICFs and from our 
Cloudy models. The solar Kr abundance is 12 + log(Kr/H) = 3.28 + 0.08 
(Asplund et al. 2005). 

^Equation (1) is used to determine ICF(Kr) for all objects except IC 418, 
where Equation (2) is used (due to the small Ar++/Ar fraction and hence 
large uncertainty in the corresponding ICF). 
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Table 12. Comparison of ICF- and Model-Derived Se Abundances 



PN 


n(Se=^+)/ 




12+ Log(Se/H) 


12+ Log(Se/H) 


Name 


n(H+)xlO^ 


ICF(Se) 


ICF 


Model 


IC 418 




22.26±8.56 






IC 5117 


2.20±0.24 


2.32±0.96 


3.71±0.20 


3.67+0.05 


NGC 6572 


0.69±0.05 


2.32±0.71 


3.20±0.14 


3.11+0.08 


NGC 6741 


0.59±0.15 


8.12±3.62 


3.68±0.25 


3.39+0.12 


NGC 6790 


0.56±0.06 


2.50±0.73 


3.15±0.14 


3.10+0.05 


NGC 6826 


0.47±0.15 


2.14±0.97 


3.00±0.27 


2.89+0.13 


NGC 6884 


1.43±0.22 


2.12±0.50 


3.48±0.13 


3.45+0.10 


NGC 6886 


0.93±0.11 


3.43±2.03 


3.51±0.30 


3.59+0.11 


NGC 7027 


1.51±0.11 


3.23±1.93 


3.69±0.30 


3.78+0.12 


NGC 7662 


1.43±0.32 


4.64±2.47 


3.82±0.29 


3.50+0.10 



Note. — Comparison of Se abundances derived from ICFs and from our 
Cloudy models. The solar Se abundance is 12 + log(Se/H) = 3.33 + 0.04 
(Asplund et al. 2005). 
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Table Al. Possible Identifications for the A5 346.0 Feature 





Wavelength 




Other Multiplet 


Ion 


(A) 


(eV) 


Lines (A) 


[Kr IV] 


5346.0 


2.318-0.0 


5867.7 




CI 


5345.6 


11.165-8.847 


5350.6, 


5356.5 


SII 


5345.7 


17.386-15.068 


5320.7, 


5346.1 


Fell 


5345.7 


12.786-10.467 


5323.7, 


5351.5 


cm 


5345.9 


42.170-39.852 


5337.4, 


5341.5, 5348.8 








5353.2, 


5359.9 


N II 


5345.9 


29.237-26.919 


5344.4, 


5347.4, 5371.5 








5374.6, 


5388.6, 5416.2 


Fell 


5346.0 


12.818-10.500 


5353.4, 


5401.1, 5407.2 








5459.6 




SII 


5346.1 


17.386-15.068 


5320.7, 


5346.1 


Felll 


5346.4 


25.417-23.099 


5314.6, 


5322.2, 5335.9 








5338.0, 


5343.7, 5354.9 








5364.9, 


5365.3, 5369.6 








5377.0 




NI 


5346.5 


14.314-11.996 


5310.4, 


5333.4 



Note. — Potential identifications of A5 346.0 are given, with 
the upper and lower energy levels of the transition and wave- 
lengths of other multiplet members. 



-54- 



Table Bl. Transition Probability and Effective Collision Strength References 



Ion Aij 

Se IV Biemont & Hansen (1987) Butler (2007)'^ 

Kr III Biemont & Hansen (1986b) Schoning (1997) 

Kr IV Biemont & Hansen (1986a) Schoning (1997) 

Kr V Biemont & Hansen (1986a) Schoning (1997) 



^In preparation. The effective collision strength of 
[Se IV] 2.287 urn at = 10^ K is 7.66. 
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Table B2. Adopted Threshold Photoionization Cross-Sections and Radiative 

Recombination Fit Parameters 





IP 


(7th 


^rr 


Element Charge 


(eV) 


(Mb) 


(10-^^ cm=^ s-^) 


Se 


9.752E+00 


l.lOOE+01 


3.120E-01 


1 


2.119E+01 


3.810E+00 


3.360E-01 


2 


3.082E+01 


1.740E+00 


4.830E-01 


3 


4.294E+01 


6.260E-01 


6.430E-01 


4 


6.832E+01 


7.870E-01 


1.220E+00 


5 


8.170E+01 


3.290E-01 


1.620E+00 


6 


1.554E+02 


2.310E+00 


2.380E+00 


7 


1.911E+02 


1.690E+00 


3.670E+00 


8 


2.270E+02 


1.260E+00 


4.640E+00 


9 


2.640E+02 


9.510E-01 


5.950E+00 


10 


3.040E+02 


7.080E-01 


7.500E+00 


11 


3.470E+02 


5.170E-01 


l.lOOE+01 


12 


3.880E+02 


3.700E-01 


1.020E+01 


13 


4.651E+02 


2.310E-01 


1.230E+01 


14 


4.750E+02 


1.510E-01 


1.490E+01 


15 


5.190E+02 


6.910E-02 


1.760E+01 


16 


6.570E+02 


3.270E-01 


1.950E+01 


17 


6.951E+02 


2.580E-01 


2.510E+01 


18 


7.360E+02 


1.950E-01 


2.760E+01 


19 


7.930E+02 


1.360E-01 


3.490E+01 


20 


8.390E+02 


8.550E-02 


3.200E+01 


21 


8.860E+02 


4.050E-02 


3.500E+01 


22 


9.870E+02 


7.260E-02 


3.940E+01 


23 


1.037E+03 


3.460E-02 


4.230E+01 


24 


2.542E+03 


1.270E-01 


4.930E+01 


25 


2.661E+03 


l.OlOE-01 


5.390E+01 


26 


2.798E+03 


7.690E-02 


5.860E+01 


27 


2.964E+03 


5.440E-02 


6.370E+01 


28 


3.123E+03 


3.440E-02 


6.890E+01 


29 


3.274E+03 


1.640E-02 


7.430E+01 
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Table B2— Continued 





IP 


O'th 




A 


Element Charge 


(eV) 


(Mb) 




/-I r\ — ^ 1 '\ —1 \ 

(10 cm"* s ^) 






3 1 1 flE- 

(J. X X L/XJ 


-02 


8 OOOF+01 


ox 


3 fi33E-l-n3 


1 48nF- 




8 5Q0F-I-01 
o . u c/ u J-j ^ X 




1 5'R7F,+n4 

X.tJ(J 1 -1— ^ 1^ VJt: 


1 4nnF- 


-09 


1 niOF+09 


33 




6.740E- 


-03 


1 n80F-l-02 




X . t:W WXJ^^ W X 


1.150EH- 


-01 


9 4PinF— 01 

^.^tJWXJ wx 


1 

X 


2 43fiE-l-ni 


5.520E4 


-00 


89nF— 01 

O.OiijWXL/ WX 


9 


0.\JCfOl-J\\J X 


5.820E-f 


-00^ 


1 OOOF-1-00 




5 24QE-l-ni 


3.850E- 


-01^ 


9 ITDF-I-OO 


4 


6 46QE+ni 


8.310E- 


-01 


1 '^fiDE+OO 


(J 


7 84QE+ni 


3.430E- 


-01 


1 690E+00 




1 iinE+n2 

X • X X v7 1 ' \ \J 


4.850E- 


-01 


9 7fiOE+00 


7 


X . ilj O O 1 Z 1 \j Zj 


2.140E- 


-01 


Q "^finE+oo 


Q 

o 


9 'ii nF,_i_n9 


1.550E+00 


A '?'?nF-l-00 


Q 


9 fiSl F-l-n9 

Z( . vjo X xi/n^ 


1.200E4 


00 


a ^ 1 DE-l-OO 


1 n 


o. xuuj— /n^uzj 


9.080E- 


-01 


7 380F-I-00 


1 1 


O . tj CJ W XJ U ilj 


7.010E- 


-01 


q n^OF-l-00 


1 9 




5.340E- 


-01 


1 inoF-i-oi 

X . xuuxjn^u X 


13 

X (J 


4 51 nE+n2 


3.970E- 


-01 


1 540E+01 


14 

Xrt 


4 Q7nF+n2 


2.890E- 


-01 


1 490E+01 

X . t: Zj V 7 1 1 \ \j X 


1 5 


5 45nF+n2 


1.970E- 


-01 


1 670E+01 


ib 


c; nQniT' 1 n9 
o.yoUrlj+U/ 


1.210E- 


-01 


i.y/uii/+ui 


17 


6.408E+02 


5.600E- 


-02 


2.230E+01 


18 


7.859E+02 


2.740E- 


-01 


2.520E+01 


19 


8.330E+02 


2.150E- 


-01 


2.890E+01 


20 


8.838E+02 


1.620E- 


-01 


3.370E+01 


21 


9.367E+02 


1.150E- 


-01 


4.330E+01 


22 


9.977E+02 


7.190E- 


-02 


4.020E+01 


23 


1.050E+03 


3.410E- 


-02 


4.210E+01 


24 


1.151E+03 


6.230E- 


-02 


4.670E+01 


25 


1.205E+03 


2.970E- 


-02 


5.030E+01 
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Table B2— Continued 





IP 








Element Charge 


(eV) 


(Mb) 




(10"" cm3 s-i) 


26 


2.928E+03 


I.IOOE- 


-01 


5.860E+01 


27 


3.069E+03 


8.760E- 


-02 


6.370E+01 


28 


3.227E+03 


6.670E- 


-02 


6.890E+01 


29 


3.381E+03 


4.770E- 


-02 


7.430E+01 


30 


3.594E+03 


2.990E- 


-02 


8.000E+01 


31 


3.760E+03 


1.430E- 


-02 


8.590E+01 


32 


3.966E+03 


2.710E- 


-02 


9.210E+01 


33 


4.109E+03 


1.310E- 


-02 


9.840E+01 


34 


1.730E+04 


1.240E- 


-02 


1.150E+02 


35 


1.794E+04 


6.000E- 


-03 


1.230E+02 



Note. — Ionization potentials, threshold photoionization cross- 
sections (Tth in megabarns (10~^^ cm^), and RR rate coefficients 
are listed. These are used in Equations (B2) and (B3) to 
determine the photoionization cross-sections and radiative rate 
coefficients as functions of photon energy and electron tempera- 
ture, respectively. 

^(7th of Kr III from equation (Bl) multiphed by 2. See §3.1. 

•^(Jth of Kr IV from equation (Bl) multiplied by 0.25. See §3.1. 
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Table B3. Adopted High- Temperature Dielectronic Recombination and CoUisional 

Ionization Fit Parameters 



Element Charge 


(cm"' s"-"- 


Ki.5) 


To 
(K) 


^di 


(K) 


P 


^ci 

(cm'^ 


-1) 


X 




K 




Kr 


3.612E- 


-09 


1.806E+04 


1.317E-10 


2.948E+05 





7.761E- 


■10 


1.169E- 


-05 


4.433E- 


-01 


1 


9.360E- 


-09 


2.894E+04 


1.129E-01 


1.052E+04 





3.260E- 


■10 


9.500E- 


-06 


4.599E- 


-01 


2 


1.555E- 


-08 


3.926E+04 


1.626E-01 


3.345E+04 





1.370E- 


■10 


1.121E- 


-05 


4.623E- 


-01 


3 


1.798E- 


-08 


4.497E+04 


2.608E-01 


8.352E+04 





5.983E- 


■11 


1.142E- 


-05 


4.567E- 


-01 


4 


1.560E- 


-08 


4.589E+04 


4.777E-01 


1.483E+05 





2.326E- 


11 


1.186E- 


-05 


4.356E- 


-01 


5 


7.783E- 


-09 


4.048E+04 


1.412E+00 


2.487E+05 





9.163E- 


12 


7.076E- 


-06 


3.378E- 


-01 


6 


5.494E- 


-09 


5.234E+04 


2.780E+00 


3.663E+05 





5.930E- 


12 


9.025E- 


-06 


3.121E- 


-01 


7 


2.018E- 


-09 


1.198E+06 


4.793E+01 


1.556E+02 





4.807E- 


■12 


8.671E- 


-06 


3.031E- 


-01 


8 


1.237E- 


-07 


1.603E+06 


4.243E-01 


1.020E+06 





5.516E- 


■12 


1.399E- 


-05 


3.015E- 


-01 


9 


1.517E- 


-07 


2.079E+06 


7.856E-01 


7.537E+05 





4.363E- 


12 


1.443E- 


-05 


3.049E- 


-01 


10 


1.787E- 


-07 


2.631E+06 


1.106E+00 


4.985E+05 





3.521E- 


12 


1.520E- 


-05 


3.083E- 


-01 


11 


8.725E- 


-09 


3.000E+06 


5.660E+01 


2.038E+05 





2.814E- 


12 


1.529E- 


-05 


3.126E- 


-01 


12 


3.803E- 


-07 


3.330E+06 


6.238E-01 


7.157E+05 





2.228E- 


■12 


1.518E- 


-05 


3.167E- 


-01 


13 


4.721E- 


-07 


3.382E+06 


5.695E-01 


1.543E+06 





1.740E- 


■12 


1.527E- 


-05 


3.187E- 


-01 


14 


5.656E- 


-07 


3.594E+06 


5.291E-01 


2.281E+06 





1.343E- 


■12 


1.517E- 


-05 


3.197E- 


-01 


15 


6.639E- 


-07 


3.736E+06 


5.028E-01 


3.156E+06 





1.031E- 


■12 


1.530E- 


-05 


3.197E- 


-01 


16 


7.642E- 


-07 


3.800E+06 


4.774E-01 


4.370E+06 





7.748E- 


13 


1.534E- 


-05 


3.178E- 


-01 


17 


1.129E- 


-06 


4.767E+06 


9.321E-01 


7.151E+06 





5.714E- 


-13 


1.556E- 


-05 


3.158E- 


-01 


18 


9.965E- 


-07 


5.182E+06 


9.855E-01 


8.913E+06 





4.381E- 


■13 


1.782E- 


-05 


3.102E- 


-01 


19 


9.745E- 


-07 


5.398E+06 


1.009E+00 


9.284E+06 





3.428E- 


■13 


1.823E- 


-05 


3.089E- 


-01 


20 


7.857E- 


-07 


4.516E+06 


1.733E+00 


1.016E+07 





2.571E- 


■13 


1.862E- 


-05 


3.051E- 


-01 


21 


4.350E- 


-07 


3.825E+06 


1.855E+00 


9.269E+06 





1.822E- 


■13 


1.895E- 


-05 


2.982E- 


-01 


22 


4.326E- 


-07 


3.600E+06 


1.091E+00 


1.222E+07 





1.285E- 


13 


1.970E- 


-05 


2.969E- 


-01 


23 


2.104E+12 


5.262E+05 


4.781E+02 


7.499E+07 





8.082E- 


14 


1.998E- 


-05 


2.914E- 


-01 


24 


5.756E- 


f04 


1.597E+03 


2.247E+07 


4.667E+05 





4.274E- 


14 


2.086E- 


-05 


2.820E- 


-01 


25 


1.018E- 


fOO 


8.355E+06 


3.344E+00 


8.891E+06 





2.159E- 


14 


1.567E- 


-05 


1.961E- 


-01 


26 


3.039E- 


-02 


1.112E+06 


1.375E+02 


1.146E+07 





1.901E- 


■14 


3.304E- 


-05 


1.618E- 


-01 


27 


4.183E- 


-02 


5.550E+05 


8.983E+01 


1.156E+07 





1.474E- 


■14 


3.370E- 


-05 


1.576E- 


-01 


28 


8.498E- 


-02 


4.947E+05 
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Note. — Fit parameters for high-temperature DR and CI, used in equations {B6) and {B7), respectively, are shown. Fits were 
performed to the DR data of Fournier et al. (2000) and the CI data of Loch et al. (2002). 
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Fig. 1. — Correlations between Kr and Se ionic fractions and those of commonly detected 
elements are plotted. Each dot corresponds to a single model from our grid of Cloudy models. 
The default PN abundances of Cloudy are assumed, and dust grains are not included. In 
some of the panels, families of curves are seen; these generally correspond to single values of 
Teff , with each curve produced by variations in nu and U. Fits to these correlations (§4) are 
overplotted as thick solid lines. 
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Fig. 2. — Same as Figure 1, except the default PN abundances of Cloudy are divided by 
ten for all elements heavier than He, and NLTE stellar atmospheres with halo abundances 
are used. This shows that our Se and Kr ICFs are minimally affected by overall (lower) 
metallicity, with the exception of the correlations shown in the two bottom panels. The fits 
shown are the same as in Figure 1. 
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Fig. 3. — Same as Figure 1, except C-rich dust is included in the models. The fits shown 
are the same as in Figure 1. The inclusion of dust physics produces negligible effects on the 
ionic fraction correlations. The same conclusion was found for PNe with 0-rich dust (not 
plotted). 



